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INTRODUCTION 


SoME years ago, the problem of hysteresis in sorption appeared to be inex- 
plicable. Recent investigations on this phenomenon in this laboratory, 
have conclusively shown the phenomenon to be real and reproducible.® *¢ 
Moreover, the interesting phenomena of drift of the hysteresis loop* and 
disappearance of the hysteresis loop? have been discovered. It has also 
been shown that the cavity concept is the general explanation of the hyste- 
resis effect and the various phenomena connected with it. In the light of 
the cavity concept the phenomenon of the disappearance of the hysteresis 
effect initially exhibited by the elastic gels has been attributed to the pro- 
perty of the swelling of such gels in solvating liquids. The importance of 
the elasticity of gels which swell on the imbibition of solvating liquids has 
been shown by studies in sorption on organo gels such as rice (a monocot),®* 7 
dhal (a dicot)? and gum arabic (a plant exudate).® The generality of the 
phenomenon has been established by further work on the sorption of water 
on certain proteins such as casein, egg albumin and gelatin presented in this 
paper. 
EXPERIMENTAL 

The McBain-Bakr quartz fibre spring technique! described in an earlier 

paper* was used in the present investigations. 


Proteins 
Casein.—Kahlbaum’s ‘‘ casein nach Hammersten ” was employed. 


Denatured casein—Hammersten casein was treated with boiling abso- 
lute alcohol for one hour and dried in an air oven at 75° C. for one and half 
hours. The product obtained was found to be insoluble in aqueous sodium 
acetate (the original casein was soluble) but soluble in alkali. 
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Egg albumin.—Metrck’s “soluble egg albumin ” was employed. 


Denatured egg albumin—A solution of egg albumin (about 10°) jn 
water was heated for one hour on a boiling water-bath, to denature the 
protein. The precipitated albumin was filtered, washed with water and 
dried at 40° C. for 2 hours in vacuum. 


Gelatin lumps.—A 20% solution of gelatin (Merck’s gold label) in boiling 
water was allowed to set in a test-tube. The gel was removed by plunging 
the tube in boiling water for a few seconds and suddenly inverting it. The 
cylindrical piece of gelatin was cut into lumps (approximately a cm. cube), 
These were dried slowly at room temperature in partial vacuum (about 
10cm. pressure) for 4 days. 


Dehydration of proteins 


In the case of casein and egg albumin, particles passing through 20- 
mesh and retained by 30-mesh sieves were used. When particles of this 
size were not available (e.g., Kahlbaum’s sample of Hammersten casein), 
the fine powder was made into a thick paste, dried in vacuum for two hours 
at the temperature of dehydration. The mass obtained was then ground 
and particles of the above-mentioned size selected. 


Dehydration of each protein was effected by heating to an appropriate 
temperature in vacuum for six hours. The temperature of dehydration was 
60° C. in the case of normal and denatured casein and 40° C. in the case of 
normal and denatured egg albumin. 


In the case of the gelatin lump, the dehydration was effected at the 
temperature of the thermostat (30° C.). 


Air-free water 


Distilled water was kept in the bulb attached to the sorption tube and 
the air was flushed out in vacuum for half an hour. 


Air-free ethyl alcohol 


Ethyl alcohol was obtained by distilling absolute alcohol twice over 
metallic calcium in an all-glass distillation apparatus. The bulb was almost 
filled with alcohol and the air in the bulb was flushed out by allowing a portion 
of alcohol to evaporate in vacuum. During sorption and desorption mea- 
surements, the possibility of the grease on the stopcock being removed by 
the solvent action of alcohol was avoided in the same way as described for 
carbon tetrachloride.® 
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Sorption and desorption of water vapour 


A suitable amount of the dehydrated protein was placed in the bucket 
attached to the spring and was degassed for 5 hours at 10-? mm. attained 
by means of Cenco Hyvac pump. After degassing, a series of sorptions 
and desorptions of the vapour of liquid on the protein at 30°C. were con- 
ducted. In all the experiments sufficient time was allowed for the attain- 
ment of equilibrium. With water, about 10 hours were found to be necessary 
for securing equilibrium and with ethyl alcohol only one hour. 


RESULTS 
Sorption of water 
Casein-water.—Casein takes at the end of the Ist, 2nd and 3rd sorptions, 
33-0 gm. 33-6 gm. and 37-4 gm. of water respectively per 100 gm. at satura- 
tion pressure. The initial hysteresis loop is reduced in size in the 2nd cycle 
of sorption and desorption and entirely disappears in the 3rd (Fig. 1). 
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Fic. 1. Casein-water (A Ist cycle, B 2nd cycle, C 3rd cycle) 


Denatured casein-water—The sorptive capacity of denatured casein 
for water, at the end of the three successive sorptions are 32-7 gm., 34:0 gm. 
and 34-3gm. of water respectively per 100 gm. The hysteresis loop 
obtained in the first cycle disappears in the subsequent cycles, Fig. 2. 


Egg albumin-water—Egg albumin takes at saturation pressure, 62-8 gm., 
68-8 gm. and 81-2 gm. of water per 100 gm. at the end of the Ist, 2nd 
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Fic. 2. Denatured casein-water (A Ist cycle, B 2nd cycle, C 3rd cycle) 


and 3rd sorptions respectively. The initial hysteresis loop disappears in 
the 2nd and 3rd cycles (Fig. 3). 
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Denatured egg albumin-water.—The sorptive capacities in the Ist, 2nd and 
3nd cycles are 38-4 gm. 45-6 gm. and 47-4 gm. per 100 gm. of the denatured 
protein. There is no hysteresis loop at all even in the Ist cycle (Fig. 4). 
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Gelatin lump-water—Sorptive capacities at saturation pressure at the 
end of the Ist and 2nd sorptions are 48-7 gm. and 55-4gm. per 100 gm. 
There is no hysteresis loop even in the 1st cycle (Fig. 5). 


Sorption of ethyl alcohol 


The sorption at saturation pressure of ethyl alcohol was measured on 
casein, denatured casein, albumin, denatured albumin and gelatin. Casein 
and denatured casein sorbed 1-8c.c. and 9-9c.c. of alcohol per 100 gm. 
of protein. Albumin and denatured albumin show no sorption at all. The 
gelatin lump takes no alcohol. 


DISCUSSION 
Swelling property and disappearance of the hysteresis effect 


The relation between the swelling property of adsorbents and the dis- 
appearance of the hysteresis effect in sorption was indicated in the studies 
in the sorption of water on rice, dhal? and gum arabic.® The present studies 
on casein, egg albumin and gelatin give additional support to this relation. 
That cavities with constricted necks entrap water during desorption and 
cause hysteresis effect has been established. In the light of the cavity con- 
cept, rigid porous adsorbents should show permanent and reproducible 
hysteresis effect on successive sorptions and desorptions and instances of 
remarkable permanence and reproducibility have already been recorded. ¢ 
With elastic gels which swell on the imbibition of solvating liquids, such 
permanence and reproducibility are not to be expected. Proteins in general 
swell in water and become elastic. Casein and albumin (normal as well as 
denatured) when dehydrated have a comparatively rigid structure. Their 
stable cavities entrap water and cause hysteresis. Near the saturation 
pressure, the gels will be in swollen condition and the cavity walls elastic. 
At this stage conditions are favourable for the disappearance of the cavities. 
The disappearance of the cavities however, is a process taking time. With 
some elastic gels the hysteresis loop disappears after the Ist cycle of sorption 
and desorpiion and with other gels at the end of the 2nd or 3rd cycle. It is 
to be noted that during desorption, fresh cavities are not formed but cavities 
still existing will remain intact. If the swelling gel is kept at saturation 
pressure for a sufficiently long time (days and even months), it is probable 
that the number of cycles of sorptions and desorptions at the end of which 
the hysteresis effect disappears is diminished. As a rule in all elastic adsorb- 
ents which swell on the imbibition of solvating liquids, there should be either 
no hysteresis loop or the loop initially exhibited must disappear on successive 
sorptions and desorptions. With nonsolvating liquids however there should 
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be permanent and reproducible hysteresis effect. This is illustrated in rice- 
carbon tetrachloride system.’ 


The observations on the gelatin lump are of great interest. In this 
adsorption experiment with a single solid lump of gelatin, free from pores, 
‘no hysteresis effect was noticed in the very first cycle of sorption and de- 
sorption. This is to be expected on the basis of the cavity concept. De- 
natured albumin behaves similarly, showing no hysteresis loop in the first 
cycle though it was activated at the same temperature (40° C.) as the normal 
albumin. Either cavities have all collapsed during denaturation or fresh 
cavities are not produced during dehydration, of the denatured protein. 


Experiments on the sorption of ethyl alcohol on casein, albumin and 
gelatin are interesting. The sorptive capacity is negligibly small. Alcohol 
is a nonsolvating liquid for casein, albumin and gelatin and there is no 
swelling in this liquid. Consequently sorption is very low with casein and 
there is no sorption at all with egg albumin. 


Another feature is the progressive increase in the sorptive capacity for 
water of the proteins on successive sorptions and desorptions. Gelatin, 
normal and denatured casein and albumin exhibit this property. Such a 
behaviour has been noticed in rice, dhal and gum arabic and is attributed 
to the loosening of the gel structure. 


DENATURATION AND SORPTIVE CAPACITY 


The sorptive capacities of denatured casein and denatured albumin are 
lower than those of normal casein and albumin. 


The subject of denaturation of proteins has attracted considerable atten- 
tion during recent years. It is a well-known fact that proteins undergo 
change in their properties on denaturation. They become less soluble and 
their swelling gets less. It is probable that on denaturation, casein and egg 
albumin become less hydrophilic and swell in water to a smaller extent and 
consequently have lower sorptive capacities than normal casein and egg 
albumin. 


The high sorptive capacity of denatured casein for ethyl alcohol (9-9 c.c. 
per 100 gm.) in comparison with the low value for normal casein (1-8 c.c. 
per 100 gm.) also indicates that on denaturation the proteins become less 
hydrophilic and show a tendency to adsorb more of the nonsolvating liquids. 


SUMMARY 


A series of sorptions and desorptions of water vapour at 30°C., on 
gelatin, casein, egg albumin, denatured casein, denatured egg albumin have 
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been conducted and in these systems either there is no hysteresis loop at 
all or the loop initially exhibited disappears on successive sorptions and 
desorptions. 


These results indicate that the swelling of the adsorbent on imbibing 
the solvating liquid is responsible for the disappearance of the hysteresis 
effect. 


Experiments on the sorption of ethyl alcohol at saturation pressure on 
the above proteins at 30°C. show that the nonsolvating liquid is either not 
adsorbed at all or taken up to a very small extent. 


With casein and egg albumin, the sorptive capacities for water are lower 
in the denatured forms of the proteins. These results indicate a decrease 
in hydrophilic character of the proteins on denaturation. 
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INTRODUCTION 


Tae problem of retention of moisture by soils has been extensively studied 
by numerous workers. Recent investigations on the phenomenon of 
Hysteresis in sorption® have thrown much light on the porous structure 
of rigid and elastic adsorbents. With a view to investigate the porous 
structure of soils, a study of hysteresis was undertaken with reference to 
the sorption of water on a red laterite soil and black cotton soil and on the 
clay, silt and fine sand fractions of the soils. The decrease in sorptive capa- 
city of the soils on the removal of organic matter has been studied. 


EXPERIMENTAL 


Mechanical analysis of soils—According to the International Society 
of Soil Science, soil passing through a 2 mm. sieve called “ Fine earth ”’ is 
divided into the following four fractions—coarse sard (2-0-2 mm.), fine 
sand (0:2-0:02 mm.), silt (0:02-0-002 mm.) and clay (less than 0-002 mm.). 
The international method A was employed for resolving the soil into these 
fractions.® In place of ammonia, sodium hydroxide was used as the dis- 
persing agent. 

Estimation of organic matter in the soil—For the estimation of organic 
matter in the soil, the new modified volumetric method using chromic 
acid®’® was adopted. The results of the mechanical analysis of the two 
soils are given in Table I. ‘The accuracy of the method is 2%. 


Activation of the soil and soil fractions.—Soil particles passing through 
20-mesh sieve but retained by 30-mesh sieve were used. In the case of soil 
fractions, the dry powder passing through 100-mesh sieve was used. The 
samples were activated by heating to 80° C. for 6 hours in vacuum obtained 
by the Cenco hyvac pump. 
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TABLE | 
| Fractions of red soil Fractions of black cotton 
Meehanical fractions per 100 gm. soil soil per 100 gm. soil 
| | 
gm, gm. 

Coarse sand “ 40-6 5-2 
Fine sand 22-6 45-5 
Silt 2-9 15-5 
Clay 30°7 33-9 
Organie matter 0-25 0-63 





Sorption and desorption of water vapour—The quartz fibre spring 
technique® was employed in the study of sorptions and desorptions of 
water vapour on the soils and soil fractions. A suitable amount of the 
activated adsorbent was placed in the bucket. The adsorbent was degassed 
by evacuation for 6 hours at 10-*mm. A series of sorptions and desorptions 
of water vapour at 30°C. was conducted. The results are indicated in 
Figs. 1 and 2. 


The sorptive capacities (for water vapour at saturation pressure) of 
red laterite soil and black cotton soil before and after treatment with 
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Sorption and desorption of water vapour on 
A. Black cotton soil with organic matter (3rd cycle). 
B. » without organic matter (2nd cycle). 
Cc. Red ‘soil with organic matter (3rd cycie). 
D. - without organic matter (2nd cycle). 
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Sorption and desorption of water vapour on 


Clay from black cotton soil (2nd cycle). 
Clay from red soil (2nd cycle). 
Silt from black cotton soil (2nd cycle). 
Silt from red soil (2nd cycle). 


hydrogen peroxide and of their fractions—clay, silt and fine sand—are given 
in Table II. 
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TABLE II 
ee 2 sal a Red soil. Sorption of water | Black cotton soil. Sorption 
Soil and soil fractions per 100 gm, of water per 100 gm, 
| 
. Nee See eee, 
gm. gm. 
Soil with organic matter a ie 8-0 18-0 
Soil without organic matter oa 6-2 | 9-9 
Fine sand oe| 0-4 0-5 
Silt i 11-2 15-8 
5 | 40-6 


Clay . 21- 








From the percentage of the various fractions of the soil and their sorp- 
tive capacities for water at saturation pressure, total sorption contributed 
by the various fractions in 100 gm. of soil has been calculated. This in 
relation to the sorption of the original and hydrogen peroxide treated soils 
has been indicated in Table III. 




















TABLE [II 
Sorption of | Sorpti f 
Red soil and its fractions P Black cotton soil and its fractions | * — ° 
water water 
Red soil original eo 8-0 % Black cotton soil original 18-0 % 
Red soil treated with hydrogen Black cotton soil treated with 9-9 % 
peroxide os 6-2 % hydrogen peroxide | 
Fine sand 7} Fine sand ] 

(22-6 gm. at 0-4%) .. 0-09 | 71-02% (45-5 gm. at 0-5%) .. 0+23 i | 16-44% 
Silt (2-9 gm. at 11-2%) .. 0°33 | Silt (15-5 gm. at 15-8%) 2-45 | | oar 
Clay (30+7 gm. at 21-56%) 6-6 J Clay (33-9 gm. at 40-6%) 13-76 J | 

DISCUSSION 


The sorptive capacities of red laterite soil and black cotton soil at the 
saturation pressure of water are 8 gm. and 18 gm. respectively per 100 gm. 
of activated soil. The difference in sorptive capacities is obviously due to 
several factors such as difference in chemical composition, soil fractions, 
organic matter and soil structure. 


As is to be expected, removal of organic matter of the soil by treatment 
with hydrogen peroxide reduces the sorptive capacity of the soil for water. 
The red soil had 0:25% organic matter while the black cotton soil had 
0-63°%. On removal of organic matter, red soil suffered a decrease in soIp- 
tive capacity from 8-0% to 6-2% and black cotton soil from 18-0% to 99%. 
The effect of a given amount of organic matter on pore space seems to depend 
on the nature of the constituents of the soil. The exact picture of the spatial 
distribution of the organic matter in soils is obscure. 
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The sorptive capacity of fine sand is negligibly small, while that of silt 
is much higher. Clay has markedly higher sorptive capacity than silt. The 
increase in sorptive capacity from fine sand to silt and from silt to clay is 
more marked in the black cotton soil than in the red soil. It is noteworthy 
that clay from black cotton soil takes up nearly twice the amount of water 
‘sorbed by clay from the red soil. Similarly, silt from the black cotton soil 
takes up more water than that from red soil, Table I] and Fig. 2. 


Knowing the percentage of the three soil fractions—fine sand, silt and 
clay and their sorptive capacities it is possible to calculate the total sorption 
of water as contributed by the three fractions, Table III. In red soil, total 
sorption of water due to these fractions in 100gm. of soil is 7-02 gm. 
whereas those of the original and hydrogen peroxide treated soils are 8-0% 
and 6:2% respectively. The total value contributed by the component 
fractions is intermediate between those of the original ard hydrogen peroxide 
treated soils. In the case of the black cotton soil, the total sorption contri- 
buted by the fractions is 16-44gm. This value is also intermediate between 
sorption by the original soil (18-0%) and hydrogen peroxide treated soil 
(9:9%). 

In all the systems investigated, permanent and reproducible hysteresis 
loops have been obtained. This indicates that unlike elastic adsorbents 
like rice,1 gum arabic,? and proteins,’ the soils and their fractions are rigid 
porous systems, having cavities with constricted ends which are responsible 
for the hysteresis effect. The soils thus retain more of water during desorp- 
tion than during sorption. 


Soil in nature is periodically exposed to extremes of weather—very dry 
and highly humid. By virtue of the rigidity of its structure, a soil subjected 
to extremes of weather, when changed from highly humid condition to one 
of normal humidity retains greater amount of water than what it would hold 
when changed from very dry weather. This excess of water amounting to 
about 194 may be negligible when soil is saturated with water in highly humid 
weather but is appreciable in weather of normal humidity in which water 
holding capacity of the soil is less than half. For example at 50% humidity, 
ted soil holds only 2% water. Compared with this, the excess 0-5% water- 
entrapped in the cavities is quite considerable. It is probable that the 
hysteresis effect exhibited by soils is of practical importance in the conserva- 
tion of moisture by soils in nature. 


SUMMARY 


Studies in sorption and desorption of water vapour at 30°C. on red 
laterite soil and black cotton soil and their fractions—fine sand, silt and 











234 Miss Sharda Gulvady and others 





clay have been carried out. Sorptive capacity of fine sand is found to be 
very small. That of silt is higher. Clay fractions have the highest sorptive 
capacity. Removal of organic matter lessens the sorptive capacity of the 
soil. 


The soils and their fractions have all yielded permanent and reproducible 
hysteresis loops. There is thus evidence that in soils we are dealing with 
rigid systems having cavities with constricted ends. The practical impor- 
tance of the hysteresis effect in soils in nature is indicated. 
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1. INTRODUCTION 


Tue problem of the dispersion of ultrasonic velocity in liquids has been in- 
vestigated by Boyle and G. B. Taylor, R. W. Wood, A. L. Loomis and 
J. C. Hubbard,” S. Parthasarathy,* Bar,* and K. G. Krishnan® All the 
investigators excepting Bar worked in the range of frequencies between a 
few Kilocycles and 20 Megacycles per secord. They observed no definite evi- 
dence of dispersion in any of the liquids studied. Bar used exceptionally 
high frequencies and measured velocities in water at 50 M.C. He too 
observed no dispersion. 


In this laboratory, it has been possible to construct a V.H.F. osefllator 
capable of reaching 50 M.C. Refined methods of comparing ultrasonic 
velocities in a liquid at different frequencies are used in a fresh attack on 
the problem. The present paper gives the results of experiments on ultra- 
sonic velocities in two liquids, namely, water and xylol in the ranges of 
frequencies 42 M.C. to 2 M.C. and 32 M.C. to 2 M.C. respectively. 


2. EXPERIMENTAL 


Water and xylol, purified by careful redistillation, are taken for investi- 
gation. The source of ultrasonics is a tourmaline plate cut in this laboratory 
and excited into resonant vibration at some one of its overtones by a 
Hartley oscillator. A specially designed Very High Frequency oscillator 
is used to work in the range of frequencies between 10 M.C. and 40 M.C. 
The circuit is the Hartley shunt-fed one with a slight modification. Filament 
chokes are used to isolate the circuit from external loops of wire. Con- 
necting leads, which at lower frequencies would serve merely to join tank 
condenser to tank coil, have greater inductance at such high frequencies 
than the tank coil itself and hence they are completely avoided. Regarding 
the tank circuit, to develop the necessary back voltage across the tank, and 
thence the high anode swing, a very low tank condenser of value -00005 is 
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used. A loop of copper strip 8 cm. in diameter, bridging the tank condenser 
serves as the tank inductance at 42 M.C. The frequency of the oscillator 
is measured using a reaction receiver in conjunction with a heterodyne wave- 
meter both of which have ranges of frequencies extending from a few K.C. 
to 15 M.C. At very high frequencies, using the method of autodyne recep. 
tion with the reaction receiver, audible beat notes are obtained at positions 
of the receiver tuning, which are submultiples of oscillator frequency. The 
receiver being kept at one of the above places, the frequency meter is coupled 
to it, and a known harmonic of the frequency meter is then made to beat 
with the oscillator. The oscillator frequency is determined from the fre- 
quency meter reading, and a knowledge of the number of its harmonic. 


The diffraction pattern is photographed on sensitive plates giving short 
exposures such that there is no appreciable temperature rise during the 
exposure. The aperture of the camera lens is cut down to less than one 
centimeter diameter to obtain an extremely sharp diffraction pattern. The 
fringe width is measured with a Hilger comparator correct to 1 to 2 parts 
in 10,000. 


The velocities are thus determined at a number of different frequencies, 
the temperature in all cases being maintained the same. Results are given 
below from which it can be seen that there is no evidence of dispersion. The 
chief source of error in the comparison of velocities at different frequencies 
is due to the slight and unavoidable temperature difference between the two 
velocity measurements. In order to overcome this criticism, measurements 
with water are made in which ultrasonic waves of two different frequencies 
are produced simultaneously in the liquid by exciting the crystal at two of 
its harmonics with the V.H.F. oscillator and the low frequency oscillator 
at the same time.* The sample pictures of the diffraction pattern with 
simultaneous excitation are reproduced in Plate XIII and it shows the 
diffraction lines due to the combination and difference frequencies n, +n, 
and n, — mg and nm, — 2ng, where n, and n, are the oscillator frequencies. The 
results reveal no dispersion of ultrasonic velocity. 


It may be noted here that at 42 M.C., the wavelength in water is so small 
that the ultrasonic grating is able to resolve clearly the full radiation from 
the mercury arc into its various components with a considerable dispersion. 
Such photographs are obtained but not reproduced here. 


* The application of the method of simultaneous excitation to the problem of dispersion 
of ultrasonic velocity was first made by S. Parthasarathy and others, Jour. Sci. and Ind. 
Research, 1945, 3, 299. 
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Ultrasonic Veloctty in Liquids at High Frequenctes 


TABLE I 


Ultrasonic velocities in water 














Temperature °C Frequency ‘Meaiaicbis Vi locity Meters 
: per second per second 
34:7 | 41-45 1540 
34°7 16-91 1542 
34°7 13-16 1540 
84-7 | 5-639 1544 
TABLE II 


Ultrasonic velocities in xylol 
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Temperature °C. 


Ultrasonic velocities in water: 
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per second | per second 
| 
32-3 31-91 | 1311 
32-3 | 20-67 1311 
32-3 16-90 1312 
32-3 9-418 1313 
32°3 2-060 1314 
TABLE III 


Simultaneous excitation 
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| per second per second 
33-8 | 39-504 | 1536 
| 9+444 1536 
33-4 39-42 | 1534 
| 13-12 | 1534 
| ! 
4. SUMMARY 


Two liquids namely, water and xylol, are investigated for the ultrasonic 
velocity in the range from 2 M.C. to 42 M.C. and 2 M.C to 32 M.C. respec- 
No evidence of dispersion is observed. 
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1. INTRODUCTION AND EXISTING EXPERIMENTAL METHODS 


Tue elastic behaviour of quasi-isotropic aggregates was treated by Voigt 
in 1887. Bruggeman (1930) gave a theory, in which an attempt was made 
to treat the boundary conditions rigorously. Special aggregates consisting 
of two kinds of isotropic particles have also been treated by Bruggeman 
(1937). These theories were tested with natural crystals such as quartz, 
fluorite, calcite, and barite; and with aggregates of these crystals occurring 
in nature. Calculated values are always higher than the observed ones. 
Birch (1938) says, “‘In any particular aggregate the question of directional 
properties may be investigated.” From this statement we can infer that he 
expects a directional variation in the elastic behaviour, but the methods in 
vogue are not powerful enough to tackle the problem. 


Adams and Coker (1906), Adams and Gibson (1931), Adams and 
Williamson (1923), Bridgman (1924, 1925), Birch (1936), Von Karman (1911), 
Roess (1935) and Zisman (1933) studied the elastic behaviour of rocks by 
static methods. Birch and his collaborators (1937, 1938, 1940, 1943), 
Ide (1933, 1935, 1937), Grimme (1938), and Kausakabe (1904-05) investi- 
gated the elastic behaviour by dynamic methods. The static methods used 
in those investigations are the common ones such as bending, twisting and 
compressing under static loads, bending and _ twisting periodically by 
impact or otherwise, bouncing spheres from a polished surface, etc. Dynamic 
methods fall into two groups according as the measurement is made (1) by 
detecting the frequency at which resonance occurs when the specimen is 
maintained in vibration by a piezo-electric or electromagnetic driving device 
or by electro-static traction or (2) by direct measurement of the frequency of 
vibration set up by tapping or otherwise exciting free vibration. All the 
above methods have certain limitations. For example, static methods are 
known to give rise to complex creep effects or elastic hysteresis. The pre- 
sence of minute cracks or other bedding planes which are either inherent in 
the rocks or produced in the preparation of the samples make the static 
method unreliable. Ide (1937) has shown that static results are always 
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lower than the dynamic ones by more than 10 per cent. “‘ Due to the inherent 
jnadequacy of the static measurements.’ Piezo-electric driving device 
requires much power and involves the complication of coupled circuits. 
The accuracy of the driving device is limited by the loading effect of the 
attached crystal. 


2. DESCRIPTION OF A NEW DYNAMIC METHOD 


Bhagavantam and Bhimasenachar (1944) have recently developed a 
new method which is of importance in studying the elasticity of all solids. 
This method, known as the “ Wedge method ”’, is based on the principle 
that a piezo-electric wedge can be made.to give out an ultrasonic beam of 
a continuously varying frequency. This method opens up new possibilities 
of geophysical investigations and for the successful application of the same 
jn determining the directional dependence of the propagation of elastic waves. 


It is well known that plane parallel plates of piezo-electric crystals when 
used in oscillatory circuits act as generators of specified frequencies, the 
actual value of the frequency being dependent on the thickness of the plate. 
Instead, when a wedge is used, the behaviour is somewhat different on account 
of its varying thickness. As the exciting frequency of the electrical circuit 
is varied, an appropriate portion of the wedge responds to the electrical 
frequency and generates an ultrasonic beam of that frequency. Thus the 
continuously varying thickness of the wedge enables a continuous band of 
ultrasonic beam being obtained. The range of the band depends upon the 
dimensions of the wedge. The lower limit is fixed by the maximum thick- 
ness and the upper limit by the power of the driving circuit, the sensitiveness 
of the detecting device and the breaking stress of the crystal employed. The 
range of the tourmaline wedge used in these investigations is from 1 to 12 
megacycles per second. For avoiding mounting effects, the area of the 
specimen under investigation must be larger than that of the wedge. This 
wedge is placed on the rock plate of known orientation. As the frequency 
of the electrical circuit is varied, ultrasonic energy which is of different fre- 
quencies when different regions of the wedge respond, passes through the 
plate. When the frequency of the wedge in a particular position coincides 
with a natural frequency or an overtone thereof of the rock section, the 
latter begins to resonate. 


The rock section and the wedge are mounted in a special holder consist- 
ing of an annular brass plate electrode secured by a brass rod to an ebonite 
plate. The rock plate of uniform thickness is placed on this electrode. The 
wedge is laid on the specimen with a film of oil in between to facilitate good 
acoustical contact. O. the top of the wedge a small brass disc is placed 
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and the whole arrangement is held pressed lightly against the brass plate 
by a fine spring. The spring acts as a second electrode. The entire set up 
is mounted on another ebonite plate and is provided with spring loaded 
levelling screws to facilitate accurate alignment. The arrangement js 
lowered into a rectang'lar cell containing carbon tetrachloride such that only 
the lower surface of the rock section comes into contact with the liquid. The 
ultrasonic beam from the wedge passes through the rock section and enters 
the liquid. A beam of light with the usual optical arrangement is allowed 
to pass through the ultrasonic grating set up in the liquid and the Debye- 

ears diffraction patterns are observed. The intensity of these patterns is 
a maximum when the rock section resonates. This permits the resonance 
position to be detected. The frequency in that position is then measured 
with an accurate wavemeter. On account of the neatness of the method 
and its decided advantage in studying even small fragments of the materials, 
the work is extended in the present studies to the investigation of the direc- 
tional variation of the elastic behaviour of representative specimens of 
sedimentary and igneous rocks. 


3. DESCRIPTION OF THE SPECIMENS 


A brief description of various rocks studied in these investigations is 
given here. 

(a) Sedimentary rocks 

Kurnool limestone-—The colour is steel black to bluish black. It con- 
tains disseminated specks of pyrites. It is cryptocrystalline. 

Bhima limestone——This specimen from the middle division of Bhima 
formation of Hyderabad state is cryptocrystalline. 

Jubbalpore marble.—It is a white medium-grained marble having a 
pearly lustre. The individual grains of calcite show interlocking texture. 

Grey Hyderabad marble.—This Pakhal marble, having a vitreous lustre, 
contains calcite, dolomite and tremolite. 

Pink Rajaputana marble.—This fine-grained pink marble has a pearly 
lustre. Calcite, dolomite and quartz are the constituents. 
(b) Igneous rocks 

Gneiss —(Uxbridge, U.S.A.). The crystals of hornblende are linearly 
arranged. Interwoven with these bands there are elongated prisms of 
zoisite. Quartz is present in abundance. The gneissic texture is clear. 

Pink Hyderabad granite-——The colour is mottled pink, white and black 


due to felspars (microperthite and plagioclase) and biotite respectively. The 
texture is hypidiomorphic inequigranular. 
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Granite of Kondavidu.—It occurs as distinct entities with defined direc- 
tional trends in the Kondavidu hill ranges. It is medium-grained holo- 
crystalline rock consisting mainly of flesh red felspar and colourless blebs of 
quattz. Microperthite, oligoclase, muscovite and magnetite are the other 
constituents. The texture is hypidiomorphic and it is inequigranular. 


Leptynite (from Kondavidu).—This rock is dominantly felspathic in 
type and is an important variety occurring along the entire range in associa- 
tion with charnockites and often in contact with kondalites. It is an equi- 
granular medium grained rock having a mottled appearance due to colour- 
less quartz blebs, pink garnets, felspars (microperthite and plagioclase) and 
biotite specks. The texture is xenomorphic granular. 


Keratophyre (Pitshead, Snowden, North Wales).—It is a porpbyritic 
rock with rectangular plates of white felspar with dark mafites represented 
by leucoxene and chlorite in a megascopically dense light coloured ground 
mass of white, bluish and greyish colour. The texture is porphyritic with 
felspathic ground mass. 


Norite of Kondapalli.—\t occurs as dykes and segregated bands. It is 
xenomorphic granular rock. Plagioclase predominates. Corcentrated por- 
tions of hypersthene contain minute flakes of brass yellow mineral (Chalco- 
pyrite ?). Diopside, magnetite and apatite are the other constituents. 


Chloritised dolerite (from Hyderabad).—It is a green rock essentially 
composed of felspar, chlorite and augite. It is very much altered. Pyrite 
occurs in cubes. The augite and the altered plagioclase are intergrown in an 
ophitic manner. 


Dunite-—{Mt. Dunn, Newzealand). It is light yellowish green on fresh 
surfaces and shows a greasy to glassy lustre. The texture is xenomorphic 
granular with predominantly olivine and chromite. 


Deccan trap from Gulberga.—This rock is dark in colour, black to greyish 
black on fresh fractured surfaces. The texture of the ground mass is sub- 
ophitic with lath-shaped felspars and interstitial augite. The texture of 
the rock as a whole can be described as holocrystalline porphyritic with 
phenocrysts of plagioclase. Labradorite and Enstatite augite are the essen- 
tial constituents with apatite and iron ores as accessories. Brown plagonite 
glass is observed. 


Basalt from Rutlam.—lIt is fine-grained (anamesite), The texture is 
glomeroporphyritic texture of Judd in which there is normal intersertal texture 
in the ground mass which contains groups of plagioclase crystals. Labra- 
dorite, augite, iron ores and glass are the constituents. 
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Eclogite (Glenelg district, Southerlands, Scotland).—These rocks jllus- 
trate the principle of minimum volume, the influence of which can be seen 
not only in the essential minerals, garnet and omphosite but also in the 
characteristic accessory mineral rutile. The retrograde change is amphi- 
bolisation. The amphibole is hornblende. The texture is granoblastic, 
Quartz is also present. 


4. RESULTS 


Plates perpendicular to three mutually orthogonal directions, taking 
easy breakage plane as reference, are cut for each rock under investigation. 
Velocities of longitudinal ultrasonic waves with frequencies ranging from 1 
to 12 megacycles along the thickness directions are determined. The 
differences of the extreme values of the velocities of each rock are expressed 
as percentages under the column “ Anisotropy” after rounding off to the 
first significant figure. Effective elastic constants are calculated and they 
are represented in the tables along with densities. There are small varia- 
tions in the densities in the same specimen for plates of different orientations. 
The variations can be attributed to the fact that the cross-sectional character 
of the units is not the same in all directions. The results on the sedimentary 
rocks are given in Table I while those on the igneous rocks are given in 
Table II. 

5. DISCUSSION OF RESULTS 


(a) Anisotropy.—Results given in Tables I and II show that there is 
an intimate relationship between the propagation of elastic waves and the 
structural trends of the rocks. The amount of anisotropy depends on the 
state of stress enjoyed by the material. 


Bhima limestone is a good example of a specimen subjected to load or 
static metamorphism. So the velocities in the two sections perpendicular 
to easy breakage directions are equal. Grey Hyderabad marble and Pink 
Rajaputana marble are impure calcareous sediments metamorphosed by 
stress and temperature. Consequently the velocities are not equal in the 
three directions. Jeffreys (1929) has suggested the possibility of sedimentary 
rocks being anisotropic due to the stratification of the layers. 


In the case of igneous rocks, the velocities in the three directions are 
different. Anisotropy increases with the influence of stress. Eclogite 
shows the least difference which is almost negligible. This is in conformity 
with the nature of its genesis namely, that it is a plutonic rock formed essen- 
tially by hydrostatic pressure at great depths. That the velosity of propa- 
gation of sound waves is different in different directions suggests that deep 
in earth’s crust where there are large differential stresses, it will be incorrect 
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to represent the elastic behaviour of rocks as that of a single homogeneous 
isotropic solid. 


As a result of extensive seismic surveys over an oil producing limestone 
in south-western Persia, two velocities were obtained by seismic methods. 
Richards (1933) explained this fact by saying that for a limestone there is no 
uniformity in elastic behaviour. It has, however, to be presumed that no 
further work was done in this direction as to a querry by Andrews in the dis- 
cussion, he replied that “‘ No evidence of the accuracy of anisotropy was 
sought as the experiments were mainly concerned with the velocities in direc- 
tions parallel to the bedding plane and plates were cut to fulfil this condition.” 


Ide’s results show that for dunite there is a change of velocity of 2% 
observed in two directions even at 4000 atmospheres. This change is attri- 
buted by Birch to porosity and he contends that at greater pressures the rocks 
behave isotropically. There is a marked variation of the velocity with the 
application of pressure. The magnitude of the observed change suggests 
change of substance not contemplated by the theory of elasticity. Porosity 
is inherent in every rock type and depends on the condition of formation 
of the rock. If this be subjected to prolonged pressures and temperatures 
of the magnitude used by the previous workers it will change its minera- 
logical composition though the chemical constitution may not change. Hence 
under such circumstances an entirely new rock may result. Thus in these 
experiments though basaltic rock behaves anisotropically, eclogite which 
has the same basaltic composition but different mineral assemblage due to 
the influence of static pressure shows the least anisotropy. In view of these 
observations, the theory that rocks behave in the same way even at low 
pressures as they do under 12000 atmospheres is no longer tenable. 


From a scrutiny of the published data with regard to the velocity of 
propagation of sound in rocks, a great variation is apparent. For example 
the reported longitudinal velocity in granites varies from 5 km.sec. to 
6-24 km.sec. One of the factors contributing to this variation is that the 
velocities calculated from the compressibility data are incomplete. This 
does not, however, explain completely the reported variation. It has to be 
recognised that the elastic properties show variation in rocks depending 
upon the directional disposition of the specimen selected. It is interesting 
to note that this range of variation is within the extreme values found for 
granites in the present investigation—4-32 km.sec. to 6-62 km.sec. depend- 
ing on the direction of selection of the specimen. Similarly the variation 
for other rocks also can be attributed to the selection of the specimen 
in different directions. 
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(b) Interpretation of seismic results——In the interpretation of seismic 
results many discrepancies are met with. Jeffreys (1929) suggest that the 
earth’s crust may be regarded as made up of thin layers of different mecha- 
nical properties with the hope that they may offer an explanation of the 
seismic observation. 

In a recent paper Gutenberg (1945) attempted to study the variation of 
the velocity of longitudinal seismic waves with the depth. He assumed 
that the velocity at a depth ‘h’ is given by 

Vi =V. + bh, 

where bis a constant varying from layer to layer. This assumption how- 
ever, does not explain the low velocities in the uppermost layers of 1 to 2 
kilometers. Gutenberg suggests that “Improvements would consist in a 
more accurate determination of the factor b and in assuming a better 
functional relationship between V and 4&4.” The results of the present 
investigation indicate that the clue lies in the second alternative. It can be 
seen that directional relationships would manifest themselves very strongly 
in the uppermost layers, thus accounting for the low velocities. 


(c) Validity of postulation of new layers.—Szveral Geophysicists have 
postulated different crustal layers from the study of the velocity of propa- 
gation of seismic waves. In view of the proved anisotropic behaviour of 
rocks, the mere differences in the velocities of propagation of seismic waves 
cannot be regarded as a sufficient criterion for postulating different crustal 
layers. 


6. SUMMARY AND CONCLUSION 


The advantages of the wedge method developed by Bhagavantam and 
Bhimasenachar for the determination of elastic properties of rock sections 
are discussed. The experimental technique underlying the above method 
is briefly described. From a systematic investigation of different rock sec- 
tions, it is shown that there is a definite anisotropy with regard to their elastic 
behaviour. In the case of sedimentary rocks, the velocity is higher in a 
direction perpendicular to the bedding or cleavage plane. The value of 
the velocity is greatest ina direction normal to the rift plane in the case of 
granites and to the easy breakage direction, in other igneous rocks. It is 
least in a direction normal to the hard way plane of granites. 


The amount of anisotropy appears to be a function of the stress experi- 
enced by the rock. The special behaviour of eclogite, in that the velocity 
is almost the same in all the directions, is in conformity with the nature of 
its genesis, namely, that it is a plutonic rock formed essentially by the hydro- 


static pressure at great depths. 
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data alone. 


It is pointed out that correct inferences cannot be drawn regarding the 
nature of the crustal layers from divergent results obtained through seismic 
The anisotropic nature of the rock through which the waves 
are propagated has also to be taken into account. The importance of the 
realisation of this fact in the interpretation of seismic data is indicated, 
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In conclusion, I have pleasure in acknowledging my thanks to Professor 
C. Mahadevan and Professor S. Bhagavantam for suggesting the problem 
and for their direction of the research. To Dr. J. Bhimasenachar I am 
deeply indebted for his constant help and advice. 


REFERENCES 


Carnegie Inst. Wash. Pub., 1906, No. 46, 69. 
Gerl. Beitr., 1931, Bd. 15, 241-50. 


Papers from Geophy. Lab. Carnegie Inst. Wash., 193), 
No. 770. 


Jour. Franklin Inst., 1923, 195, 475-579. 
Proc. Ind. Acad. Sci., 1944, 20, 298. 


Jour. Geol., 1938, 46, 59-87, 113-41. 

Bull. Geol. Soc. America, 1936, 47, 1235. 

Jour. Applied Physics, 1937, 8%, 129. 

Jour. Geol. 1940, 148, 752-66. 

Bull. Geol., Soc. America, 1943, 54, 63-86. 

Proc. Amer. Jour. Sci., 5th series, 1924, 7, 81-102. 
Amer. Jour. Sci., 5th ser., 1925, 58, 483-98. 


Elastizitat Konstanten von Kristallaggregc.ten, The Hague, 
1930. 


Ann. Der. Phys., 1937, 29, 160. 

Phil. Mag., 1935, 20, 304, 

Amer, Jour. Sci., 1945, 243 A, 285-312. 

Proc. Nat. Acad. Sci, 1933, 22, 81-91, 482-96. 
Jour. Geol., 1937, 45, 689-716. 

Rev. Sci. Inst., 1935, 6, 295-98. 


The Earth, Its Origin, History and Its Physical Consti- 
tution, 1929, 2nd ed., Cambridge, pp. 85-100. 


Zeit. D. Ver deutischen Ingenieure, 1911, 55, 1749-57. 
Jour. Coll. Sci. Tokyo., 1904-05, 20, Art. 19, 

Proc. Phy. Soc. London, 1933, 45, 70-81, 

Comptes Rendus, 1935, 200, 147-49. 

Proc. Nat. Acad. Sci., 1933, 19, 653. 





he 


Or 


im 


31, 


jue, 


Sti 











INTERACTION OF NITROUS OXIDE, NITROGEN 
AND PHOSPHORUS UNDER SILENT ELECTRIC 
DISCHARGE 


By G. S. DESHMUKH, M.Sc., AND Y. D. KANE, M.Sc. 


(Asstt. Prof. of Chemistry, Benares Hindu University) 


Received November 26, 1946 
(Communicated by Prof. S. S. Joshi, D.sc., F.a.sc.) 


Tue thermochemical behaviour of a number of gases towards phosphorus, 
particularly its glow, has been studied by various workers. There is, how- 
ever, but little quantitative information in the literature on the production 
of these changes under well-defined conditions of electrical discharge. 
Joshi? ** gave an important lead in this field of research in observing the 
general utility of data for the time-variation of the discharge current, in the 
elucidation of the mechanism of the corresponding reaction. This is illus- 
trated by numerous results in these laboratories, especially on the reduction 
with hydrogen of gaseous nitrous oxide and in frozen films at liquid air 
temperature.»* It was of interest therefore, to extend these observations 
to a study of the reaction between the phosphorus films and nitrous oxide, 
under silent discharge. 
EXPERIMENTAL 


The general arrangement of the apparatus and circuit employed are 
shown in Fig. 1. A Siemens’ type ozoniser made with soft glass tubes 
(fitted co-axially with a rubber greased ground-joint) was surrounded by a 
tube filled with salt solution, to serve as a low tension electrode. The whole 
system was immersed in a jacket, in which a stream of water was kept flowing 
at a constant rate in order to minimise the heat effects under the discharge. 
The ozoniser was connected through tap 1 to a mercury manometer and 
through tap 2 to a Tépler. Nitrous oxide obtained from a commercial 
cylinder was admitted through tap 8 and stored in glass reservoirs R, and R, 
which had been evacuated previously. With taps 8, 10 and 5 closed, on 
opening tap 7, most of the nitrous oxide was frozen to a solid of negligible 
vapour pressure in traps cooled by liquid air in Dewar vessels D, and D,. 
The uncondensible gases were removed through tap 5 by the Topler, tap 2 
being closed. Dewar vessels D; and D, were then carefully removed; the 
first portion of the frozen gas to vaporise was pumped out; tap 5 closed; 
and the middle fraction stored in R, and Rg after toplering off through tap 
5 the last part of the condensed gas to vaporise. After the entire assembly 
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of the apparatus was tested for vacuum for at least 48 hours, taps 3 and 1 
were closed. The inner electrode of the ozoniser was then coated uniformly 
on the outside with a concentrated solution of phosphorus in carbon di- 
sulphide, and fitted in the ozoniser system. The solvent was removed first 
with a hyvac pump from tap 3 and finally on the Topler. At the commence- 
ment of every experiment nitrous oxide from R, and R, was frozen in traps 
cooled by liquid air in D, and D,; the middle fraction of the vaporised gas 
was led over a train of tubes charged with P.O, and KOH, and introduced 
into the annular space of the ozoniser through tap 3 at the desired pressure. 


The A.C. obtained by means of a 1 K.W. Rotary Converter worked off 
220 volt, D.C. mains was fed to the primary of a 3 KVA transformer; one 
of its secondaries was earthed; the other connected serially through a 10,000 
ohm stabilising resistance, was dipped in the salt solution forming the inner 
electrode of the ozoniser, The outer electrode was earthed through an 
oxide rectifier type A.C. micro-ammeter uA (see Fig. 1). The applied 
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potential expressed in kilo-volts (r.m.s.) and called V, was calculated from 
a knowledge of the stepping up ratio of the transformer and of the primary 
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potential, V,; the latter was regulated by means of a variable resistance 
jn its primary. 
DISCUSSION OF RESULTS 


Series of experiments started at various gas pressures were made on 
the nitrous oxide-phosphorus reaction under conditions mentioncd above; 
a fresh film of phosphorus was employed in each series. In agreement with 
a general finding due to Joshi,” ** on changes under electrical discharges 
especially in gases, the interaction between nitrous oxide and phosphorus 
sets in only above a characteristic ‘ threshold potential ’ V,,, where the material 
in the annular space breaks down as a dielectric. V,, depends, inter alia, 
on the gas pressure p and frequency of the electrical supply. The course of 
the reaction was followed by observations at regular intervals of the gas 
pressure p at constant volume, and also of the corresponding discharge 
current i. These data for one typical series correspording to a nitrous oxide 
pressure of 12:2cm. Hg, and 4 kilo-volts (r.m.s.) for the secondary potential 
of 50 cycles frequency are shown in Fig. 2. In all cases observed, p decreased 
progressively during, and was negligibly small at the erd of, the reaction. 
A further continuation of the discharge even at enharced V the secordary 
potential had no effect as judged from the constancy of the spectral nature 
of the glow and of the gas pressure p. As shown by the total time for the 
reaction and especially by the average time-rate of the pressure decrease, 
the rate of the change increased with V and 1/p. This follows readily from 
Joshi’s general finding that the ‘threshold potential’ V,, increases with p 
and that the rate of a discharge reaction and the correspording i deperd 
(increase with) on V-V,,.% 4 


As observed by Joshi” ** in the decomposition of (pure) nitrous oxide 
under silent discharge, the current i varied sensibly discontinuously during 
the present reaction at a given V, indicative of consecutive changes (see Fig. 2). 
Thus, for example, during the short initial stage when p diminished, i ircreased 
appreciably. In the next stage which was the longest, both these quantities 
varied comparatively slowly. The last stage is characterised by a rematkably 
sudden and rapid decrease of i to a constant minimum at which p also becomes 
negligibly small, suggestive of the completion of the reaction. From results 
of a number of reactions at various initial p and V, of which the curves in 
Fig. 2 are one typical illustration, it was found that the last stage during 
which i dropped suddenly, commenced when the residual p of the reaction 
mixture reached 3 to 6mm., i.e., about 4 per cent. of the initial p. This 
sudden drop of i at a constant V, was also observed when nitrogen was used 
instead of nitrous oxide (Fig. 3). 
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The observed decrease in the over-all p (as also the corresponding 
changes) are explicable, on the following mechanism for the reaction: 
N,O>N,+0O (ia) (The atomic products in \ia), (ib) and (ii) are 
N,O > NO +N (ib) removed immediately by phosphorus). 


NO >~N+0O0 (ii) 


Changes (ia) and (ib) occur simultaneously and represent the initial 


stages ; (ia) is suggested to be the more predominant reaction. 


with Joshi’s results for the decomposition of (pure) nitrous oxide under 
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silient discharge ;** mixed with an excess of hydrogen;> and when frozen 
over glass and certain metallic receptacles at the liquid air temperature.* 
Atomic oxygen and nitrogen produced in (ia) and (ib), react with the pbos- 
phorus film, the former practically instantaneously. The fact that free oxygen 
was not detected in the reaction mixture at any stage is in accord with this 
deduction. So far as the gas phase is concerned, the net change represented 
by this stage (ia and ib) is 2N,O - NO + Ng since the other (atomic) pro- 
ducts, viz., N and O being removed by phosphorus do not contribute to the 
gas pressure. It follows from this, that the gas pressure during this stage of 
the reaction should be constant. This anticipation was realised experi- 
mentally by working on the following suggestion due to Prof. Joshi:* It is 
to reduce the over-all velocity to a minimum by maintaining just V,, on the 
ozoniser, so that the otherwise rapid and even transient reactions are slowed 
sufficiently to enable observations of their time-rate and also the nature of 
the corresponding products, which ordinarily would escape detection. With 
this restraint a definite though short lived stage during which p remained | 
sensibly constant was observed initially in a number of N,O — P reactions 
under the discharge. It may be emphasized here that the discontinuous 
time-variation of i characteristic of the stage-wise progress of this reaction 
(Figs. 2 and 3) was more conspicuous when V did not exceed V,, appreciably 
as is to be expected from Joshi’s views.* Studies of other types of discharge 
reactions in progtess in these laboratories have shown the wide applicability 
of this device, viz., excitation at V,, for elucidating the characteristic course 
and mechanism of a discharge reaction.® 


Since according to Joshi,” ** V-V,, determines i, its rise in the initial 
stage (see Fig. 2) may be attributed to a reduction of the corresponding V,,,, 
due firstly, to N,O changing to NO and nitrogen; and secondly, due to 
phosphorus being replaced by its oxide (and possibly nitride, vide infra), 
This deduction agrees with the observation by various workers that the oxida- 
tion of phosphorus produces ionisation in the surrounding air;? and also 
with our experimental results with phosphorus films reacting with air (dried 
carefully with the arrangement shown in Fig. 1) subjected to the ozoniser 
discharge at 27-2, 31-4, 47-9, 53-4, 60-0, 65-7 and 73-2cm. Hg pressures 
at V varied from 1 to 5kV of 50 cycles frequency, and the time of exposure 
to discharge being increased up to 3-5 hours. It was always found that 
whilst p decreased, the corresponding i increased markedly. 


The decomposition of nitric oxide according to (ii) constitutes chiefly 
the next stage. Its atomic products would be removed from the gas 
phase, oxygen especially being affected very much faster due to its much 
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greater affinity towards phosphorus. The production of nitric oxide was 
detected in a series of experiments at various initial pressures on nitrous 
oxide (a fresh film of phosphorus being used in every case) by introducing 
from tap 3, a small bubble of air in the ozoniser after exposure to discharge 
for about 3-5 minutes; this reaction mixture when led over a small amount 
of pure mercury (shown by Hg, in Fig. 1) fouled it immediately due to for- 
mation of nitrogen dioxide; this was never observed after longer exposures; 
and also in absence of air when the mixture was toplered out at any stage 
of the reaction. Nitric oxide occurs therefore, in the present reaction as 
an intermediate product during only early stages. This conclusion was 
confirmed further by analysis of mixtures from reactions at various initial 
pressures of nitrous oxide in the range 15 to 42cm. Hg; the discharge was 
interrupted after different periods of exposure. The mixture was led by 
opening tap 1 and lowering the manometric mercury in M, through a trap 
cooled by liquid air in D, in order to separate nitrous oxide; nitric oxide 
also freezes to a solid with a small sublimation pressure and can therefore, 
be toplered out completely from tap 2. The mixture thus pumped out con- 
sisting of nitric oxide and nitrogen was collected in a gas burette and analysed 
for nitric oxide by absorption with acidified ferrous sulphate solution. In 
no case was nitric oxide observed after 10 minutes i.e., after about 1/30th 
of the total time for the reaction; furthermore, the ratio NO/Ng in the mixture 
increased as the time of exposure was reduced. Nitrous oxide does not 
react with phosphorous thermochemically “even when touched with a piece 
of red hot iron.”*® At higher temperatures it inflames, “‘ the products being 
phosphorus pentoxide and pure nitrogen dioxide.”’* This last is also observed 
when nitric oxide is substituted for nitrous oxide in the thermal change. It 
is not unlikely that nitric oxide is formed as an intermediate product at very 
carly stages of the high temperature reaction between phosphorus and nitrous 
oxide by analogy with the present results for the discharge reaction. The 
latter differs from the purely thermal change, in that nitrogen dioxide is not 
observed at any stage. Presumably, this applies also to the interaction under 
discharge, with nitric oxide (vide infra). 


A series of experiments was next carried out in order to investigate the 
removal by phosphorus of nitrogen from the gas phase, produced by the 
decomposition of N,O and NO according to (ib) and (ii). The ozoniser 
in Fig. 1 was filled with cylinder nitrogen (after bubbling slowly through 
traps containing freshly prepared alkaline pyrogallol in order to remove 
traces of oxygen) at 7:0, 8-0, 9-0, 11-3, 12-7 and 21-5cm. pressures ard sub- 
jected to discharge at V varied in range 2-4kV at 50 cycles frequency for 
periods up to 12 hours. In every case p diminished and the correspondin g 
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j increased. Towards the end of the reaction, i dropped remarkably 
suddenly; curves in Fig. 3 refer to a typical run with 11-3cm. of nitrogen 
at4kV. The general resemblance of the progress of this reaction, as shown 
by the time-variation of i and p with that, when nitrous oxide (initial p and 
applied V being smaller) was used instead, is also striking. For the nitrous 
oxide-phosphorus and nitrogen-phosphorus reactions at the same potential 
viz. 4kV and similar initial pressures 12-2cm. N,O and 11-3cm. Ng (see 
Figs. 2 and 3), the total time for the reaction, i.e, when the gas pressure 
became negligibly small, was 340 and 330 minutes respectively ; this similarity 
issuggestive. As with nitrous oxide, after its final sudden drop to a minimum, 
on increasing V, i first increased (almost) momentarily to a large value; this 
remained stationary as long as the corresponding V was constant (see Fig. 3). 


If after the removal of nitrogen by phosphorus, a further emount of it 
was introduced, p and 1/i continued to decrease, but on the whole more 
slowly than before. Thus e.g., in one experiment at 4kV and initial p of 
11-4cm. Hg of nitrogen, 50 per cent. of the pressure. change occurred in 
160 minutes; the total time required was about 330 minutes. At the end 
of this period, when p became negligibly small (see Fig. 3) fresh nitrogen 
was introduced at 4cm. pressure in the ozoniser. It now took about 290 
minutes for the same, i.e., 50 per cent. pressure change, despite a larger re- 
action rate due to 1/p and V (6 kV) being greater than before. With 21-5cm. 
pressure due to a further addition of nitrogen, exposures to discharge at 
10-7 kV for 39 hours produced only 60 per cent. of the total pressure decrease. 
Detailed experiments showed that any thermo-chemical change between the 
phosphorus film and the nitrogen (also air) in the ozoniser was negligibly 
small compared with that produced under the discharge. It was also noticed 
incidentally during this stage, that the ‘threshold potential’ V,, character- 
istic of the aerial oxidation of phosphorus by an ozoniser discharge is 
defined more sharply, if the inner electrode supporting the phosphorus film 
is cooled (to about 3° C.) than at the room temperature viz., 27°C. It was 
also found that changes in the appearances of the glow from initial violet 
to yellow; and of the electrode film of phosphorus from white to yellow, 
brown and finally almost black, were comparatively more distinct. 


Results for the nitrogen-phosphorus reaction under discharge due to 
fields at 500 cycles and also much higher frequencies obtained with a serial 
spark, were found to be on the whole similar to those observed at 50 cycles 
frequency (Fig. 3). 


Examination under microscope of the products (of a number of re- 
actions) on the inner electrode after long exposures to the discharge showed 











254 G. S. Deshmukh and Y. D. Kane 


small crystals of (unused) yellow phosphorus mixed with a brown substance, 
The mixture on treatment with water gave ammonia when treated with 
Nessler’s reagent, indicating that a phosphorus nitride was produced under 
the discharge. Using Geissler type discharges and metallic electrodes, 
Strutt® (now Lord Rayleigh) found that phosphorus absorbed nitrogen, and 
that the gas was not reliberated on heating, which he attributed to the forma- 
tion of a chemical compound; this agrees with the results of Newman,® 
and Moldenhauer and Dérsom.” In the present work (vide supra), the 
reaction products always revealed a nitride; its formation would appear 
therefore, to represent by far the main change in the P—N, interaction under 
the silent discharge. Judged by the time-rate of the variation of both 
p and corresponding i, this ‘clean up’ of nitrogen by phosphorus was more 
rapid as the applied V and 1/p were increased; ultimately, whole of the 
input nitrogen was removed. As observed in the last stage of each of the 
reactions with nitrous oxide also (see Fig. 2), when the pressure of the residual 
nitrogen was very small, i decreased remarkably suddenly. This supports 
the deduction made that the decomposition under discharge of nitric oxide 
is completed during early stages; and that the ‘clean up’ of nitrogen is 
practically the only change, which constitutes a later and by far the longest 
stage in the N,O—P interaction under discharge. The remarkably sudden 
drop of the current i (cf. Figs. 2 and 3) to stationary values at the end of 
P—N,O and P—N, reactions may be attributed to the production of a vacuum 
condition in the discharge space by the final ‘clean up’ of nitrogen by the 
phosphorus film, which is practically the same sole change during the latter 
part of both the above reactions; incidentally, this illustrates the marked 
‘getter’ action of phosphorus. The current i is constant as long as V is con- 
stant e.g., at 4kV in Fig. 3; it, however, increases immediately to a higher 
constant value (cf. the step-wise, current-time curve, continued from A to A’ 
in Fig. 3) owing to increased ionisation of the residual gas (and phosphorus 
vapour) under fields due to larger V. 


On general considerations, it is shown easily that the discharge current 
iin a gas under a constant applied field would increase as p is reduced, since 
the energy gained by motion under the field by an ion (chiefly electron) at 
the end of a given free path increases with its length i.e., with 1/p; the limi- 
tation, due to the circumstance that the total number of electrons leading 
to ionisation and therefore, contributing to i also diminishes with p, is not 
large. At moderate pressures, obtained during the main course of the re- 
action of N,O and nitrogen with phosphorus, i is maintained as a balance 
of two factors mentioned above; on the whole, i tends to increase as the 
pressure is diminished (see Figs. 2 and 3), except when it becomes negligible 
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towards the end of both these reactions. The large though momentary 
rise of i on applying 8 kV before the attainment of the final stage (Fig. 3) in 
this nitrogen-phosphorus interaction requires, however, further investigation. 


In conclusion, we welcome this opportunity of expressing our sincere 
and grateful thanks to Professor S. S. Joshi, D.sc. (London), F.a.sc., Head of 
the Chemistry Department, for suggesting the problem, instructive advice 
and kind encouragement during this work. 


SUMMARY 


Interaction with phosphorus films of nitrous oxide, nitrogen and oxygen 
in silent discharge under fields due to potentials (V) of 50, 500 cycles and the 
high frequencies has been studied at various pressures p. In agreement with 
a general finding due to Prof. Joshi,” *** the change sets in above a minimum 
V,, and that both the current i and the reaction rate at a given V depend 
upon V-V,,. p decreases progressively to zero, when the reaction termi- 
nates; this is characterised by a remarkably sudden drop of the correspond- 
ing i. Free O, is not observed; NO occurs only at early stages. Two 
simultaneous changes (ia) and (ib), the former being the more rapid, followed 
by (ii) represent the reaction mechanism viz., N,0 ~ N,+ O(ia) N,JO— 
N+NO (ib) NO>N + O(ii). The atomic products react with phos- 
phorus immediately. Removal of molecular N2, from (ia) and that formed 
by recombination of N from (ib) and (ii) by phosphorus, is the longest stage. 
Phosphorus nitride is the chief product with both N,O and N,. The time- 
variation of i is markedly synchronous with, but more pronounced than 
that of p; and more helpful in revealing the stages of the reaction, suggested 
as a general proposition by Joshi.* Interaction with oxygen produces a 
rise of i in striking contrast with that with nitrogen. 
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1. INTRODUCTION 


As is well known, when a beam of plane polarised light traverses a crystal 
of quartz along the optic axis, the plane of polarization is rotated conti- 
nuously. The rotation increases very rapidly as the wavelength of the light 
diminishes: it is 26°, 42° and 148° per millimetre respectively for the 
wavelengths A 5461, 44358 and 42537. If the crystal is viewed by an 
observer in a direction transverse to the beam, the electric vector in the 
light beam would be alternately parallel and perpendicular to the line of 
sight for each successive rotation of 90°. Hence, if the light diffused by 
the crystal is partially or completely polarised, the track of the beam as seen 
by the observer should exhibit alternate minima and maxima of intensity, 
the distance from one dark band to the next being the distance of travel of 
the incident beam in which there is a rotation of 180°. The widths of the 
bands for the wavelengths A 5461, A 4358, and A 2537 would be 7:06, 4-33 
and 1-22 millimetres respectively. A phenomenon of this kind was 
observed and photographed in smoky quartz with visible light by the present 
Lord Rayleigh (1919); it had its origin in the Tyndall scattering by the 
inclusions present in the crystal to which the smoky colour is ascribed. The 
present paper records the observation of the phenomenon in perfectly colour- 
less quartz having its origin in the diffusion of light by the atomic vibra- 
tions in the crystal lattice. Both visible light 14358 and the ultraviolet 
radiation A 2537 of the mercury arc in quartz have been successfully used 
to photograph the effect. Owing to the very great intensity of the resonance 
radiation emited by the water-cooled, magnet-controlled mercury arc 
which was employed, as also the large scattering power arising from the 
\~* law and the high rotatory power of quartz in the vitraviolet, the 
phenomenon is much more striking as observed with A2537 than with 
A 4358 radiation. 
2. EXPERIMENTAL TECHNIQUE 


For the purpose of the experiment, a perfectly colourless and trans- 
parent sphere of quartz crystal 5-24 centimetres in diameter was employed. 
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As was shown long ago by Sir C. V. Raman (1922), such a sphere in which 
no visible inclusions are present exhibits a uniform blue track wken 
traversed by a condensed beam of sunlight. The use of a sphere is conve- 
nient as light can be sent in any required direction. To reduce the 
surface reflections, the sphere can be immersed in clean distilled water. 
When the beam of condensed sunlight traverses the sphere along its optic 
axis and the blue track is viewed transversely, the scattered light shows 
very marked polarization in which the electric vector is perpendicular to 
the optic axis. Again, when unpolarised sunlight is sent perpendicularly 
to the optic axis and the track is observed in a direction transverse both 
to the optic axis and to the beam of light, the scattered light is strongly 
polarised with most of the light having its electric vector parallel to the 
optic axis. However, when the observation is made along the optic axis 
with the light beam perpendicular to it, the scattered light is completely 
unpolarised. This is evidently due to the fact that scattered light before 
reaching the observer traverses the crystal along its optic axis and hence 
its plane of polarisation is rotated. The amount of rotation varies for the 
different wavelengths of sunlight and the emerging light hence appears to 
be completely unpolarised. 


The foregoing observations indicated that the phenomenon of the 
bands could be expected and that for photographing them monochromatic 
light would have to be used. Another important condition for the success- 
ful performance of the experiment is that light should traverse the crystal 
along its optic axis pretty accurately. To secure this, the sphere is placed 
between two polaroids, and observed in mercury arc light and adjusted 
until the melatrope is observed centrally in the field of view with rings 
around it. By rotating the analysing polaroid clockwise, the rings of any 
particular colour appear to expand and hence the sphere rotates the plane 
of polarisation anticlockwise as viewed by the observer. The sphere is 
hence a left-handed crystal. 


To photograph the phenomenon of the bands, the radiation of a water- 
cooled, magnet-controlled mercury arc was employed as the source of light 
and an aperture of 5/8” in diameter placed just in front of the arc limited 
the beam of light used. Actually, the whole radiation of the mercury arc 
was used, as the other wave-lengths were weak compared to the A 2537 
radiation and were also scattered comparatively feebly. The light beam 
was condensed by means of a 17cm. lens of 4cm. aperture into a brass cell 
fitted on opposite sides with two strain-free vitreous silica windows of 3.cm. 


aperture. The inside of the brass cell was painted black, and the quartz sphere 
A4 
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was placed inside it on an ebonite ring and was immersed in clean distilled 
water. The light before entering the brass cell was polarised by a large 
nicol which transmits the A 2537 radiation. The sphere was adjusted so 
that the light traverses it along its diameter exactly parallel to the optic axis 
and there was thus no refraction of light at the surface of the sphere. The 
light emerging vertically from the surface of water was photographed by a 
camera fitted with a 6-5cm. quartz lens of aperture 3cm. and in this way 
a large quartz window for the observation was dispensed with. The camera 
and the top of the brass cell. were wrapped up in a black cloth to prevent 
stray: light entering the camera directly. A process plate was used and an 
exposure of one hour was given.. A photograph of the track in the trans- 
parent quartz sphere showing a large number of bands was obtained, an 
enlarged reproduction of which is given at the end of the paper. Another 
photograph with a different arrangement, using the visible 4358 A line of 
the mercury arc was also obtained, but a long exposure of 100 hours was 
needed for this purpose. Even so, the photograph was not dense enough 
to be satisfactorily reproduced. 


3. REMARKS 


The width of the bands was measured and the magnification produced 
by the combination of the sphere, water and the camera was calculated. 
From‘this the actual width of each band was deduced. This distance should 
correspond to a rotation of 180° for the particular wave-length used, and 
hence enable the rotation per millimetre to be calculated for the two wave- 
lengths A 2537 and A 4358 respectively. The results are given in the table 
below along with the standard values for comparison. 


es TABLE | 





| 
Actual width | Calculated | Rotation per 
of each band | rotation per (mm. in degrees 


| \ 
| Measured average 


Wavelength | \idth of each Magnification 








in | band in mm.® | of the system in mm. mm. in degrees |from I. C. T. 
| 
2537 1-421 (22) | 1-340 1-085 | 165-9 148-6 
} 
4358 4-133 (3) 0:97 4-25 42-4 41-6 





* The number in the brackets gives the total number of bands measured whose average 
width is given. 





The fair agreement between the calculated and the actual rotations and 
the uniform clarity of the bands indicates the highly monochromatic 
nature of the light scattered and that the other radiations are not recorded. 
The contrast between the maxima and minima of intensity of the bands 
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shows that the scattering of light in quartz is well polarised. The picture 
of the bands reveals at a glance the continuous rotation of the plane of 
polarisation taking place inside quartz and the marked polarisation of 
transversely diffused light in quartz arising from the atomic vibrations in 
the crystal lattice. 


As will be noticed, the bands are not straight, but are slightly curved. 
This is because the surface at which light enters the quartz crystal is curved, 
and hence the points of equal rotation lie along a curve parallel to the 
surface. 


In conclusion I wish to thank Sir C. V. Raman for suggesting the 
problem and for his constant encouragement. I wish also to take this 
opportunity to express my grateful thanks to Dr. R. S. Krishnan for his 
kind help which was always available. 


SUMMARY 


When a plane-polarised beam of the intense A 2537 radiation from a 
water-cooled, magnet-controlled mercury arc is sent along the optic axis 
of a perfectly clear and transparent sphere of quartz crystal free from 
inclusions, the track when photographed from a direction transverse to the 
beam exhibits striking fluctuations in intensity along its length. The 
effect, which is reproduced in the paper, is obviously connected with the 
rotation of the plane of polarisation of the polarised 2537 radiation as 
it traverses the crystal. The distance between one dark band to the next 
corresponds closely to a rotation of 180° of the plane of polarisation of the 
\2537 radiation. A similar effect was observed and photographed in 
smoky quartz by the present Lord Rayleigh in 1919 using the Tyndall 
scattering of visible light by the inclusions in the crystal. In the present 
case the effect is due to the diffusion of light arising from the atomic 
vibrations in the crystal lattice, and the clearness of the bands indicates that 
such scattering is strongly polarised. 
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1. INTRODUCTION 


As is well known, quartz is a uniaxial crystal exhibiting optical activity and 
its behaviour in respect of radiations of specified wave-length to which it 
is completely transparent may be specified by three constants, viz., its 
ordinary refractive index m,,, extraordinary refractive index n,, and rotatory 
power p (along the optic axis), in degrees per millimetre. The accurate 
measurement of these quantities has been the subject of many investigations 
and various attempts have been made to express the experimental data and 
their dependence on the radiation frequency by appropriate formule. The 
experimental work and the formule proposed up to 1927 are reviewed 
in Sosman’s The Properties of Silica (1927) and Lowry’s Optical 
Activity (1935). 


Besides the literature reviewed in these treatises some contributions 
to the subject may also be mentioned here. New formule have been pro- 
posed by Coode-Adams (1928), Bradshaw and Livens (1929) and Engl 
(1936). There has also appeared a detailed analysis of the infra-red absorp- 
tion of quartz by Saxena (1940). Servant has extended the measurements 
on optical activity to 1525A.U. (1938) and determined the birefringence 
up to 1450 A.U. (1937). The spectroscopic behaviour of quartz has been 
investigated in detail by Tousey (1940) in the remote ultra-violet. 


The formule so far proposed give values which steadily deviate from 
the newer measurements with diminishing wave-length. This is perhaps 
not surprising since the absorption frequencies in the older formula 
(especially in the ultra-violet) are empirically assumed ones. Accordingly, 
in the present investigation, the author has tried to derive new formule 
which are free of these defects. 
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2. THE FORMUL2 PROPOSED 





These are 
0003952 0003808 0003114 
_ 4 Pisin 2 sci oa = a 
(1) M* = 4°57877 + \2—9-014161 * 4 —0-011236 + 1® —0-000974 
44-224 713-55 
+ y= 78-22 * \®— 430-55 
00042647 0-0041290 . 0-0028067 


—0-014161 * A® — 0-011236 * A® — 0-000974 
41-678 . 701-354 


+ S773 * WF — 388 


(2) nt = 4-730561 + 43 


pe Ae 467 23110000815 
(3) P = x39 -014161 + A®— 0-011236  A® — 0-000974 (A®— 0-0195)? 
— 01905 


(A, in all cases, is in microns.) 


These formulz are largely refinements on the work of Coode-Adams 
andLowry. The first two ultra-violet terms, with peaks at 1190 A.U. and 
1060 A.U. replace their single term at 1127 A.U. which is nearly the mean 
of these two. The justifications for their introduction lies in the work of 
Tousey, who has found (1) an absorption cut off at 1436A.U. exhibited 
by even thin slices, (2) a sharp, strong peak of selective reflection at 1190 A.U. 
and a diffuse peak about 1060 A.U., both remaining unchanged in wave- 
length and changing only slightly in intensity for different crystal orienta- 
tions. The use of these terms with nearly equal strengths in place of the 
empirical term of Lowry and Coode-Adams extends the validity of the 
formule into the remote ultra-violet. 


The third term at 310 A.U. appears in the formule of Lowry ard Coode- 
Adams and is unambiguously indicated to be necessary by .the refraction 
data. Experiments have not advanced far enough into the ultra-violet to 
verify the existence of a specific maximum absorption at such a wave-length. 
Theoretically, however, the transitions of the L-electrons of silicon (with its 
neon-like core) to vacant outer levels would give rise to a term with a 
wave-length of about 200A.U. and the transitions of the L-electrons of 
oxygen to a term above 500 A.U. The result would be to give a lumped 
effect as if there were an absorption at some intermediate wave-length. 
It will be noticed that the relative importance of the terms at A1190 and 
41060 in the formulg is approximately the same for the ordinary and extra- 
ordinary rays, but the relative importance of the term at A310 is distinctly 
less in the case of the extraordinary ray. The manner in which the 
birefringence alters with the frequency is influenced by this circumstance. 

















262 T. Radhakrishnan 





In the case of the ordinary ray the dielectric constant and the infra-red 
absorption frequencies have been taken without change from an old formula 
of Rubens which considers the lattice vibrations strongly active in respect 
of the ordinary ray (viz., one at 21, and three about 9 which have been 
lumped together into one term). A similar formula has been constructed 
by the author for the extraordinary index, taking into account the lattice 
oscillations strongly active for the extraordinary ray (namely one at 194 
and two others near about 9, which have been lumped together). The 
appropriate value of the dielecttric constant and the oscillator strengths were 
determined from the refraction data. 





The formula for optical activity was first designed to be asymptotically 
equivalent to the Lowry and Coode-Adams formula at long wave-lengths. 
The resulting large errors in the extreme ultra-violet were found to be 
removed by introducing a square term with a very small numerator, and 
with a characteristic wave-length below 1400A.U. This might be an 
effect of the continuous absorption below A 1436, whose contribution to 
refraction is negligible. The constant 0-1905 appearing in the formula 
may possibly be the effect of a far-off infra-red term whose variation 
in the range of -14y to 3, is negligible. The term at A 310 is of less im- 
importance in optical activity than in refraction asis to be expected, since py 
optical activity is more connected with the outer inter-atomic linkages than 
the inner structure. 
3. COMPARISON WITH FACTs* 
(1) The Ordinary Ray _ 


Wave-length | Refractive Index | Refractive Index 








(microns) (Calculated) | (Observed) 
0-185467 1-67589 1-67578 
0+274867 1-58765 1 -58752 
0+396848 1-55812 155813 
0-486133 1 -54965 1-54968 
0+387563 1 -54423 1 -54433 
0-670786 154135 1-54147 
0+ 794763 1 -53835 1+53848 
1-00 153489 1-53503 
2-0582 1-51981 1-51998 
3-000 1-49947 | 1-49962 
4-09 1-4615 | 1+4620 
5-0 1-410 1-417 
7-0 1-095 | 1-167 

33 2-39 2-53 
52 2+22 2+23 
63 2-19 2-19 
83 2-167 | 2-15 
co 2-14 2-11 





* Values of refractive indices and rotatory. power above 0-185 are taken respectively 
from the treatises of Sosman and Lowry quoted before, Sosman’s estimate of the more 
probable values of nm, and m, in air and at 18°C. being used. Data below 0-185 4 are 
taken from Servant’s paperts. 


= 

















Wave-length 


| 


| Refractive Index 


‘Refractive Index 


(2) The Extra-ordinary Ray 
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(microns) (Calculated) (Observed) 
| 
0-185467 | 1-69007 1-68997 
0-274867 | 1-59819 1-59813 
0-396848 1-56776 1-56772 
0-486133 | 1-55900 1+ 55898 
0-587563 | 1-55343 1-55343 
0-670736 | 1-55047 1-55048 
0-794763 | 1-54741 1-54739 
0-84467 1 -54639 1-54640 
1-0000 | 1-54381 1.54381 
1-8000 1-53246 1-52342 
2-0582 1-52814 1-52814 
2-5000 | 1-5194 1-5195 
3.0000 | 1-5070 1-5070 
33 2-40 2-63 
52 | 2-25 2-30 
63 2-23 2 +25 
83 | 2-20 2-22 
co | 2-175 2°15 
(3) Birefringence 
Wave-length Ny te He — Ny)X10* | (we — 2) X10 
(Microns) (Calculated) (Calculated ) (Calculated) {| (Observed) 
ee rg eee ee ee ee ico alia ‘a { 
0-14546 | 1-85975 1-88365 | 23-90 | 23 64 
0-15478 1-78722 1 -80709 19-87 19-57 
0-16619 1-73169 1-74863 16-94 16-72 
6-17680 | 1-69703 1+ 71224 15-21 | 15-05 








Wave-length 
(microns) 


0-1525 
0-1550 
0-1575 

0- 1625 
0-1650 
-1675 
+1700 
185398 
+ 198979 
+ 226909 
+ 232749 
+ 273955 
- 3693999 
*435834 
520908 
-579066 
-6707846 


coococoocoeooro 


-~ 


342 
- 430 
400 
3-210 


wee 


{ 
| 





767 
722 
678-5 
600 
564°5 
533-5 
504 
371-3 
295 +33 


28+ 237 
22-5425 
16-536 


3-857 
3-132 
1-073 
0-512 





(4) Optical Activity 


Calculated rotatory;Observed rotatory 
power (deg./mm.) | power (deg./mm. ) 


= i a 


7176 
724 
678-5 
599-5 
564°5 
533 
504-5 
370-9 
295-65 
200-9 
187-247 
122-123 
60-058 
41-551 
28-245 
22-5465 
16-5359 


3-894 
3-168 
1-074 
0+520 





—_——— 

























T. Radhakrishnan 


























Th 
4. DISCUSSION 
cee : C 
The formule proposed are valid with fair accuracy over the whole range 
of measurements, while the older formule, highly accurate up to the near 
ultra-violet, fail badly in the remote ultra-violet. The following tables 
(the first one taken from Servant’s paper on the birefringence) illustrate 
the deviations of the older formule. 
a l | I 
Wave-length (A.U.) | 2210 | 2099 1932 1515 | 1464 
(%e — wy) X 10° observed a 11-97 | 12-45 13-46 20-72 | 23-64 
| | 
(%¢— my) X 10% calculated | 11-96 | 12-41 13-32 18-59 | 90.39 
from Coode-Adam’s formulz| | | 
ELLE S: TP Cerne | 
<< eee ee ee 
Wave-length ..| 0 2174 d 1700 | » 1550 
Observed p | 226-91 | 604-5 | 724 
p Calculated from Lowry | 226-98 509-5 | 748 


Coode-Adams-formula | 
This progressive error is almost completely annulled in the formule 
proposed by the author and this suggests that the changes he has made 
are at least in the right direction. The wave-lengths 41190 and A 1060 are 
known correct only to 10 A.U., and the employment of more precise values, 


when these are secured by experiment, might improve the fit. 


I voice here my gratitude to Sir C. V. Raman for suggesting the 
problem and his interest in the investigation. 





5. SUMMARY 


New formule using the latest data available and based on the observed 
spectroscopic behaviour of quartz have been given for the refraction and 
optical activity of quartz, which fit the facts over a wider range of frequencies 
than any formula proposed so far. 


APPENDIX 


The older formule from which the present ones have sprung are set 
out below :— 


Rubens’ formula for the ordinary ray :— 


) 0-0010654 , 44-224 , 713-55 
Mu* = 4° S787 + ya 9.010627 + 278-22 t ¥® 430-55 
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Coode-Adams’ for the ordinary and extraordinary rays :— 


0008067 0-002682 127-2 

2 5 lenis on sitont i eee ee ee 
| my® —~ 3°53445 + yx ~9.0197493 + 2 —0-000974T 42 108 
| ng? 3-5612587 4. 0700844614 0-00276113 127-2 


A2— 0-0127493 * A® — 0-000974 © A®— 108 
Lowry and Coode-Adams’ formula for optical activity :— 
9-5639 2-3114 


°~ 8 0-0127493  r®— 0-000974 9" 1905 
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1. INTRODUCTION 


A stupy of the temperature-variation of the refractive indices of solids is 
interesting from many points of view. On the experimental side, optical 
measurements are capable of very high precision and thus it is possible to 
obtain accurate data. Theoretically, a satisfactory theory of the variation 
has not been forthcoming so far. In fact, it is found that many crystals 
like rocksalt, fluorspar and quartz exhibit a diminution of their refractive 
indices for visible light on heating, while: others like diamond (Sella, 1891) 
and zinc-blende (Maria Mell, 1923) show an increase. Also, some of the 
crystals of the former type show the reverse type of behaviour, namely an 
increase in the refractive index with temperature, in the extreme ultra-violet. 
An adequate theory of the phenomenon must be capable of explaining all 
these facts. 


The experimental methods of measurement are also capable of improve 
ment. The usual procedure in finding the variation of refractive index 
with temperature (dn/dt) is to determine the absolute index (m) at various 
temperatures and then to deduce dn/dt from these measurements. It 
would, of course, be much more accurate if one could devise a method in 
which dn/dt is directly measured as a function of temperature. Such a 
method has been designed by the author and is discussed fully in a later 
section of this paper. The method has been applied to the case of diamond 
and dn/dt has been measured over the range of temperatures from —180°C. 
to + 450°C. for three wavelengths in the visible region. 


A phenomenological theory of the temperature-variation of the refrac- 
tive index has also been developed in this paper. It is based on the idea 
that the refractive index can be expressed in terms of the number of 
dispersion centres per unit volume and the polarisabilities of these centres. 
The latter again is dependent on the dispersion frequencies of the system. 
The variation in the refractive index of the solid can thus be represented 
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by two terms, one arising from the change in the number of dispersion 
centres and another from the variation in the spectral position of the dis- 
persion frequencies. The idea is applied to the case of diamord in this 
paper and it is shown that the rate of change of the dispersion frequencies 
with temperature is correlated to the rate of change of other electronic and 
lattice frequencies observed for the crystal. The application of the formule 
to the case of fused silica enables one to calculate dn/dt over the whole 
spectral range, which are found to fit remarkably well with the experimental 
data. This will be reported in Part II of these series. 


Thevmo-Optic Behaviour of Solids—1 


2. EXPERIMENTAL PROCEDURE 


The principle of the method may be described as follows. One employs 
a plate of thickness / and obtains interference fringes by the light obtained: 
from the two surfaces of the plate at normal incidence. (Either the trans- 
mitted or the reflected system may be employed.) The path difference 
between the two interfering beams is x = 2nl, where n is the refractive 
index. Sippose one measures the change in this path, Ax for a known 
change of temperature Az. Then, 


Ax/At =2n. Al/ At +21. An/ At, 


so that 
Ani Ax_,. 1A! 
at 2 @&t 1 At 
= (1/21) (Ax/ At) — na, (1) 


whefe a is the coefficient of linear expansion at temperature ¢. Thus, if a 
is known, then by measuring Ax and At, dn/dt can be accurately obtained. 
Of course, n has also to be known, but only to 3 significant figures for 
getting an accuracy of 1° and not 5 or 6 significant figures as are necessary 
with the usual prism methods. Another advantage of the interference 
methods is that the value of dn/dt obtained is for the absolute refractive 
index and that no corrections are required for the change in the index. of 
air aS is necessary with the prism methods. 


In practice, two types of experimental arrangement can be visualised 
employing the above principle, one or the other being suitable for a parti- 
cular case. Suppose that one is able to obtain plane-parallel plates of the 
specimen, as is for instance possible with mica. It is best in this case to 
form Haidinger rings with monochromatic light with the plate kept in a 
suitable furnace. When heated, the rings expand or contract and a 
Measurement of the expansion or contraction gives one the alteration in 


path. Onae method of obtaining Ax/At is to accurately measure Ax for 
Aba 
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a known difference of temperature At. But if one is interested in measur. 
ing dn/dt over a whole range of temperatures it is much mofe convenient 
to make a continuous observation. The temperature of the furnace js 
raised steadily at a very slow rate and the temperature corresponding to 
the passage of successive fringes across a fiduciary cross-wire can be noted, 


If optically parallel plates are not available, another method can be 
used, which was in fact employed with diamond in the present experiment. 
The crystal used for the purpose was N.C. 73 from the collection of Sir C. Y. 
Raman, which was a rectangular plate with a thickness of 0-670 mm. 
The surfaces of the plate were deposited with a thin layer of aluminium by 
the evaporation process in order to increase the reflecting power. When 
viewed by transmission in a parallel beam of monochromatic light, localised 
interference fringes were observed, which consisted of sharp bright almost 
circular rings. These were produced because of the variations in the 
thickness of the specimen. A fine line on the surface of the crystal formed 
a reference mark across which the fringes shifted as the crystal was heated. 
It may be mentioned that this method of observation is the same as was 
employed by the author in determining variations of the absolute path of 
a ray of light through the crystal plate for calculating the photo-elastic 
constants of diamond (Ramachandran, 1947). 


The diamond was kept in a suitable furnace and the fringes were 
observed by means of a: microscope. The continuous observation method 
mentioned in connection with the Haidinger rings was employed with these 
fringes also. The rate of heating was so slow that it took 6 hours to heat 
the crystal from room temperature to 400°C. The temperature was 
measured by a standardised thermocouple and the readings of the measur- 
ing instrument were such that an accuracy of +1°C. could be obtained 
even at the highest temperatures. The experiments were performed with 
three different wavelengths, viz., the 4358 and 5461 A.U. radiations of the 
mercury arc and the yellow radiations of 5893 A.U. from the sodium lamp. 
In the first case, more than 40 fringes shifted across the reference mark in 
heating from room temperature to 400°C. and in the last, nearly 30. The 
observations were repeated with the cooling crystal, this also taking roughly 
6 hours. The temperature readings taken during heating and cooling 
never differed by more than one unit in the last decimel place, i.e., nearly 
1° near about 400° C. and less at lower temperatures. For each wavelength, 
the experiments were also repeated reversing the positions of the source 
of light and of the observing microscope. This was done in order to correct 
for any movement of the crystal or of the whole furnace due to heating. 
Actually, no difference could be noticed between the two readings. Six 
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independent experiments were performed with each wavelength and the 
finally calculated values were the means of these. 


The values of a at various temperatures were available from the 
measurements of R. S. Krishnan (1946). Knowing At corresponding to 
each wavelength shift, dn/dt was calculated for the mean temperature ¢. 
When these values were plotted and a smooth curve drawn through them, 
it was found that the points appreciably diverged from the mean. This was 
obviously due to the restricted accuracy of the temperature measurement. 
Consequently, the values of dn/dt were later calculated from the temperature 
difference corresponding to a shift of 3, 4 or 5 fringes. These fell very well 
on a smooth curve, as will be seen from the curves reproduced in Fig. 1. 
Here, the continuous curves represent the mean value of all the experiments 
for each wavelength. The points denote the values from a single experi- 
ment calculatea in the above manner. 








iL A. iL. A 
“ae o 100 2 200 300 
Mone. a -7-___so é Cc 


a—6 43584 =e 054614 uu 58 93A 
Fig. 1. Curves showing the variation of dn/dt with temperature. 
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The measurement of dn/dt was also made for temperatures down to 


— 180°C. Here, the accuracy was not as high as at higher temperatures, . 


since the rate of variation was small and only 4 fringes shifted over the whole 
range and also since measurement could be made only during warming and 
not during cooling. dn/dt was calculated from the temperature difference 
for the shift of one fringe at a time. The crystal was kept in a copper 
block which was screwed on to the bottom of a brass tube in which was kept 
liquid air. The tube with the crystal was enclosed in an outer jacket of glass 
and the intervening space was kept evacuated. The thermocouple was kept 
in intimate contact with the bottom of the brass tube. The whole arrange- 
ment was first cooled to the temperature of liquid air and then allowed to 
warm up. It took nearly 3 hours for the system to attain the room tempe- 
rature. During this period, temperatures were measured corresponding to 
the passage of each fringe across the reference mark. It was roticed that 
upto about — 150°C. or so there was no shift at all of the fringes, showing 
that the rate of variation of the refractive index was extremely small at these 
temperatures. 


An obvious disadvantage of the method described here is that it is not 
very convenient for work in the ultra-violet; but this may be remedied by 
the use of photography. The method is, however, particularly useful when 
one requires to find dn/dt accurately over a range of temperatures. 


3. RESULTS OF THE EXPERIMENT 


The experimental values of dn/dt at various temperatures for the three 
wavelengths are shown by the continuous curves in Fig. 1. The values 
read off from these ( x 105) are shown in Table I and will be useful for 
reference. These are probably accurate to 1%. 




















TABLE I 
rc. | 4358 | 5461 | 5893 °C. | 4358 | 5461 | 5893 
-150 | 0-02 | 9.01 | 0-01 160 2-33 2-04 | 2:00 
-140 | 0-05 | 0-03 | 0-03 180 2-60 2-26 2-22 
-120 | 0-13 | O10 | 0-09 200 2-78 2-45 2-41 
-100 | 0-24 | 0-21 “| 0-19 | 220 2-88 | 2-57 2-53 
— 80 0-38 0-32 | 0-30 240 2:96 2-66 2-62 
— 60 | 0-53 | 0-45 | 0-48 || 260 3-06 2-76 | 72 
— 40 0-68 | 0-59 | 0-57 || 280 3-20 2-88 | 2-84 
— 20 0-83 | 0-73 | 0-70 300 3-35 | 3-03 «| 2-98 
0 | 0-99 0-87 | 0-84 || 390 3-50 3-17 | 3-12 
20 1140S | si 0-98 || 340 3-65 3-30 | 3-25 
4o | 1-299 | 1615 1-12 360 3-80 3-43 | 2 3+38 
60 | 144 | 1-299 | 1-26 380 3-94 3-56 | 3-50 
80 | 1-59 1-43 | 1-40 400 4-08 3-68 | 3-62 
100 1-75 1:57 | 1-54 420 4-20 3-80 | 3-74 
120 1-91 1-71 | 1-68 440 4-31 3-91 3-84 
140 2-07 1-85 | 1-82 | 
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The most interesting fact that comes out of the investigation is that 
dnjdt is positive for diamond over the whole range of temperatures. This 
has been observed by Sella (loc. cit.) from room temperature to 90°C. His 
value of + 1-86 x 10-5 for a mean temperature of 50°C. and 5893 A.U. 
is however too large, the present experiment indicating a value of about 
41:20 x 10-5. It is found that dn/dt is extremely small for temperatures 
below — 150°C. and that it increases practically linearly with temperature 
from about — 100° to 450°C. However, a slight hump is observed in the 
graph of dn/dt vs. t in the neighbourhood of 210° C., as may be seen from 
Fig. |. Very accurate and repeated sets of readings have been taken in 
this range of temperatures and it has been verified that the hump is a genuine 
effect. The deviation of dn/dt is however only of the order of 10%. In 
fact when the values of n calculated from the data for dn/dt are plotted 
against ¢, it is practically impossible to detect the kink in the curve, even when 
plotted on a highly magnified scale. The author is not in a position to 
give a definite explanation for the hump at 210°C. It may however be 
pointed out that a similar variation, though in a much more prorourced 
manner, is shown by the refractive index of fused silica in the neighbourhood 
of 575°C., the temperature of a—f transformation of quartz. Here, the 
kink is observable in the graph of refractive index against temperature. 
Itis not clear what kind of transformation, if any, occurs in diamord near 
about 200°C. It has however been observed (Chandrasekharan, 1946) 
that diamond exhibits very bright thermoluminescence at about this tempe- 
tature. Similarly, Nayar (1941) found that the intensity of luminescence 
diminishes rapidly at these temperatures and practically vanishes at 300°C. 
The particular diamond that was studied exhibited a feeble blue luminescence. 
‘It would be interesting to repeat the dn/dt measurements with highly 
luminescent and non-luminescent diamonds to discover if the magnitude 
of the hump has any relation to luminescence. Another interesting fact 
is that a particularly low value for the Raman frequercy of diamord has 
been observed near this temperature (R. S. Krishnan, 19464). In Fig. 1 
of this reference, it will be noticed that the two values at temperatures 
213°C. and 242°C. fall well below the continuous curve through the other 
points. 


Leaving aside such minor differences, it has been calculated that the 
tefractive index of diamond 7, at temperature t can be represented by a 
formula of the type 

Nz == No (1 + at + Bt?) (2) 
from about — 100°C. to 400°C. The values of a and § for the three wave- 
lengths studied are given in Table II. 
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TABLE II 
Wavelength A, v. | a x 108 B x 108 
4358 9-9 7-8 
5461 | 8-7 7-0 
5893 8.4 7-0 


4. THBORY OF THE THBRMAL VARIATION OF REFRACTIVE INDEX 





Pockels (1902) has recognized that the rate of change of refractive 
index with temperature (dn/dt) can be divided into two parts, namely, a pure 
temperature change which we may denote by d/dt and another part due to 
the change in density which accompanies the temperature change, 
(dn/ dp), - dp/dt. Thus, 

& _(*) 2 4(2 

dt ph, at ts, 
Of these, the quantity A may be calculated from photo-elastic data. For 
a cubic crystal, 

A=—2,. Gut2pu) 7 - 
where y is the coefficient of cubic expansion and p,, and py, are the elasto- 
optic constants. The same relation is true for isotropic solids like glasses, 
while similar relations hold for crystals of other classes. When A is 
actually calculated in this manner, it is found that it is not in general equal 
to dn/dt. Pockels, in the paper cited above, has shown that A is always 
algebraically less than dn/dt so that B is a positive quantity for a number 
of glasses of varying composition. In Table LII below, dn/dt, A and B are 
given for a number of crystals. The data have been collected from the 
Handbuch der Physik, Vol. XXI, 1929, Ch. 19, except for diamond for which 
the data are taken from a paper by the author (1947). 


=A +B say (3) 











TABLE III 

| 
Substance dn/dt+10* A-10¢ B-10¢ 

| 
Rocksalt —37°4 —37-2 ~— 0-2 
Fluorspar —12°1 —13-5 + 14 
Quartz-ord — 6-50 —12-80 + 6°30 
» «= ext ~ 7°54 — 14-10 + 6°36 
Calcite-ord + 0-60 - 8-9 + 0-6 
«6 ext oe +10°57 - 19 +13°5 
Topaz-a oe + 8-1 — 0-38 + B85 
” -B ee + 9-0 - 0-713 — 9-7 
oo. = ee + 8-0 — 0-684 + 8-7 
Diamond +10°5 + 38 + 6°7 
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It will be seen that, in all the cases except rocksalt, B is positive. Even 
for rocksalt one cannot say that the negative value, which is less than 1% 
of the observed dn/dt is not spurious. We may thus say that, in general, 
crystals and glasses exhibit a pure thermal coefficient of variation of refractive 
index which is positive. 


This picture does not however go deep enough into the problem except 
for separating the density and temperature effects. In fact, one can give 
a quantitative theory of the thermo-optic behaviour by considering the 
dispersion of solids and its origin. It has been the attempt of many workers 
to fit the dispersion data of solids into a formula of one of the following 
types :— 

nt = C’ +ZA,/(a* — 2,3) — Ki? (5) 


n? = C + ZB, A*/(A? — A,*) — KA? (6) 


where the term KA? may sometimes be absent. These two forms are really 
equivalent and can be converted one into the other by merely changing the 
constants. Also, the last term arises from the presence of an infra-red 
frequency or frequencies whose positions are not known, so that in general 
we may write (6) in the form 


n® =C + Za, A? A, — *) = C +Za,v,* — v) (7) 


where v is the wave number of a radiation of wavelength A. Now, the 
Ketteler-Helmholtz theory (and also the more recent Kramers-Heisenberg 
theory) of dispersion leads to a formula of the type 


m®=1+La/(,2 — v) (8) 


where all the dispersion frequencies are taken into account. Usually, 
the empirical formula is derived from observations in the ultra-violet, visible 
and infra-red regions of the spectrum. If there is a dispersion frequency 
in the remote ultra-violet, it will contribute a constant term to the right- 
hand side of the formula, which will then be of the type (7). We shall, 
however, consider only a formula of the type (8) in our theory since it is 
of the most general type. It is obvious that the quantities a, in the numerator 
are proportional to the number of dispersion centres per unit volume and 


‘to the oscillator strength of each centre. 


Now, if such a formula is true for a particular temperature, then a 
similar formula with slightly different constants must hold for a slightly 
different temperature, so that one can obtain an expression for dn/dt by 
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differentiating these expressions with respect to t. One thus obtains from (8) 


dn * 2a,v, dv, 1 _ da, 
ne 2 [- Gace a tos EI 9) 


Since a, is proportional to p, the density of the substance, da,/dt = (a,/p), 


(dp/dt) = — ay, where y the coefficient of cubical expansion.* Thus, 
the second term in the right-hand side of (9) becomes 
Z—y:a,(v -- v*) -= — y(n? — 1) from (8). 
dn yin®—1) | a,v, — dv 
. 2n n = on —v) at (10) 


This is a general expression for dn/dt in terms of known constants of the 
substance. v, is known from the dispersion formula, but dv,/dt is not known 
in many cases. But it may be calculated from the values of dn/dt for a 
particular wavelength and then formula (10) can be verified to hold over 
the whole range of wavelengths. It will be convenient to express dv,/dt 
in the form — x, »,, whefe x, is just a constant expressing the rate of change 
of wavenumber as so many per wavenumber of the dispersion frequency. 
It may be noted that x, is also equal to (1/,) (dd,/dt). Then, 


aS _wee-h i a, v,? ] 

dm Tae GF REX | 
—_ (11) 

—_ —yint—1) 1s a, XM, | 

2n n* Q? — 432 * | 


In this expression, dn/dt has been split up into two parts, one due toa 
change in the number of dispersion centres per unit volume and another 
due to the changes in the dispersion frequencies. We may denote these 
by P and X respectively. The term X arises from the change in the dis- 
persion frequency with temperature; but it is not at all necessary that this 
change should be a pure temperature effect. It is necessary to distinguish 
between the change that is produced by the accompanying change in volume 
yor density) and the change due to temperature alone. Jn symbols, 

m ow iw, mM. ww 
dy w dt’ w dt 

=Q +R (say) (12) 
Obviously, P arises completely from the change in density, so that one gets 
immediately, 


A=P+Q B=R (13) 





* This implies the assumption that the oscillator strengths, i.e., the transition probabilities 
of the dispersion centres, do not alter with temperature. 
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These considerations are, of course, purely formal, but they are’yet interest- 
ing. It is seen that one has to take into account a possible shift of the 
dispersion frequencies even ta explain the alteration of refractive index with 
pressure. The shift of the frequencies with pressure is not negligible. In 
fact, for diamond at 25° C. and 5893 A.U., we have from (11) P= — 3- 1 x 10-* 
and combining this with the data in Table I, we have 


P = —3-1 x 10% Q= + 6-9 x 10% R = + 6:7 x 10° (14) 


Here, the shift in frequency, due to the change in volume is as much as that 
due purely to change of temperature. 


The considerations put forward above refer strictly only to the 
behaviour of cubic crystals or of isotropic bodies. In other crystals, the 
linear thermal expansion is not the same in all directions and the refractive 
indices in different directions are influenced differently. We shall not go 
into these questions in this paper. 


5. DISCUSSION OF THE RESULTS OBTAINED WITH DIAMOND 


We shall now apply the considerations put forward in the previous 
section to the case of diamond. The dispersion of diamord can be repre- 
sented by the formula 


gy AAPA a, = 3-859 x 10, ay = 1-080 x 10° - 
~~ FAA, =10600A.U., A. = 1750 AU. ( 


(Peter, 1923). Consequently, dn/dt is given by Eq.(11). It will be noticed 
that both the dispersion frequencies are in the extreme ultra-violet and that 
one of them at 1060 A.U. is nearly 40 times as strong as the other at 1750 A.U. 
We shall make the assumption that y, = xg, in other words that both of 
them vary in proportion to their frequencies. It will be seen presently that 
this assumption is not unjustified. Besides, since A, is very weak, slight 
differences in the value of yx, will not matter. Thus, putting xy,= x2=x 
its value calculated from the measurements of dn/dt for the green and yellow 
radiations comes out to be 7:6 x 10-* at 50°C. It is interesting to note 
that this is of the same order as the value of x for the fundamental lattice 
vibration of 1332cm~. From the data of R. S. Krishnan (1946 4), it has 
been calculated that y at 50°C. is 9-8 x 10-* for this frequency. It is thus 
found that the rate of variation of the ultra-violet dispersion frequencies 
is nearly the same as that of the lattice frequency in the infra-red provided 
the rate is taken as so mamy wave-numbers per wave-number per degree. 
In fact, this behaviour is exhibited by practically all the electronic lines 
found in the absorption and emission spectrum of diamond. In Table III 
the wavelengths of various such lines are given both at liquid air and room 
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temperatures and the total change is calculated and given as so many 
wavenumbers per 1000 wavenumbers (1000 A ¥/¥) in the third column. 

















TABLE IV 
X at liquid air A at room -- 
temperature temperature ane Gaye Reference 
Disperson Frequencies 0-6 This paper 

2306 0-9 ] 
i a +K. G. Ramanathan (1946) 
2360 2364 1-8 J 
3034 3037 | 1-0 } 

3067 1-0 
3154 3157 | 1-0 + Robertson, Fox and Martin (1934) 
3201 3204 1-0 
3298 3302 1-3 J 
3447 3450 | 0-83 ) 
4152 4156 0-83 
5032 5038 1-2 
5115 5123 1-6 
5198 5204 1-2 , 
5286 5292 1+] + A. Mani (1044) 
5359 5367 1°5 
5451 5455 0-7 
5695 } 5701 1-0 
5758 5768 1-7 
1333-2cm.~* 1331 -8cm.~? 1-1 R. S. Krishnan (1946 4) 








The approximate constancy of the numbers in the third column is 
significant. It shows that the behaviour of the electronic lines is similar 
to that of the lattice frequency. This is not surprising when one remembers 
that the electronic lines arise from transitions of the very electrons that are 
responsible for the binding forces in the lattice. 


The refractive indices at various temperatures were calculated from the 
data for dn/dt by numerical integration assuming the values at 20°C. given 
by Peter (loc. cit.). In Fig. 2 the course of the refractive index-temperature 
curve has been plotted from these calculations. It will be noticed that the 
refractive index does not at all change at very low temperatures. In fact 
below — 100°, the change does not even amount to one unit in the fourth 
decimel place, which is the accuracy to which absolute indices ate known. 
Extrapolating to the absolute zero and taking Peter’s value for the refractive 
indices at 20° C., the indices at absolute zero for the three wavelengths have 
been calculated to be 


2-4490 at 44358, 2-4225 at A5461, 2-4168 ai A 5893. 
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An attempt was made to calculate the course of the refractive index- 
temperature curves, assuming that the dispersion frequercies vary in the 
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Fic. 2. The variation of refractive index of diamond with temperature. The three 
curves are from top to bottom for wavelengths 4358, 5461 and 5893 A.U. respectively. 


same way as (a) the lattice frequency or (5) the electronic line at 4152 A.U. 
Data for the temperature variation of these are available from R. S. Krishnan 
(1946 b) and P. G. N. Nayar (1941) respectively. The experimental values of 
(nm, — ngs) are shown along with those calculated from assumptiors (a) and 
(6) in Table V below. It will be seen that although the general trerd of the 
curve is shown by the calculated values, the numerical agreement with 
experiment is not very satisfactory. This only shows that the variation 
With temperature of the dispersion frequency does not exactly follow that 
of the lattice frequency or the other electronic frequency at 4152A.U. 
although it is of the same order, 
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TABLE V 
Temperature ‘c.| — 180 | 25 100 | 200 300 400 
"7 l 
M-n25 (Exp.) ..| —0-0008 0 +0-0009 | +0-0026 | +0-0055 | +0-0098 
do from (a).:| —0-0021 0 +0-0015 | +0-0038 | +0-0063 | +0-0082 
do from (4)..| —0-0014 0 +0-0017 | +0-0052 | 




















It was mentioned above that x was found to be 7-6 x 10-* at 50°C. 
A recalculation of dn/dt from this value of x showed that it fitted with 
experimental values for 5893 and 5461 A.U., but was too small for 4358 A.U. 
The calculated values were 1-20, 1-22 and 1-27 x 10°*, while experiment 
gave values of 1-19, 1-22 and 1-36 x 10-* respectively for the three wave- 
lengths. The author is of the opinion that the discrepancy for 4358 A.U. 
is due to the presence of a weak absorption at 4156A.U. in the diamond. 
It may be pointed out that on account of its proximity, the contribution 
to du/dt by this frequency at 4358 A.U. is nearly 20 times those at 5461 and 
5893 A.U. It will be interesting to see whether the discrepancy is larger 
with a diamond having a stronger absorption at 4156A.U. and also to 
actually detect the presence of anomalous dispersion of refractive index 
near this wavelength. 


My sincere thanks are due to Prof. Sir C. V. Raman for the kind 
interest that he took in the investigation and for the many discussions I had 
with him. 

SUMMARY 

A phenomenological theory has been put forward for the variation ot 
the refractive indices of solids with temperature, based on the idea that the 
refractive index can be expressed in terms of the number of dispersion centres 
and the polarizabilities of these centres, which ate again dependent on the 
dispztsion frequencies. The refractive index variation can-then be repre 
sented as the sum of two terms, one arising from the change in the number 
of dispersion centres and the other from the variation of the frequencies. 
Accurate data of the thermal variation of the refractive index of diamond 
have been obtained for the first time from — 180°C. to 450°C. for three 
wavelengths, 4358, 5461 and 5893A.U. A new method was employed for 
this purpose which measures dn/dt directly and is capable of general applica- 
tion. When applied to diamond, the theoretical ideas bring out the inter- 
esting fact that the rate of change of the dispersion frequencies with tempe- 
fatufe is correlated to the rate of change of the fundamental: lattice 
frequency of the crystal. 
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1. INTRODUCTION 


In Part I of these series, some general considerations were put forward 
indicating how a phenomenological theory can be developed for the varia- 
tion of refractive indices of solid substances with temperature. In this 
paper, a quantitative application of the theory is made to the case of fused 
quartz. It is shown that the theory can account satisfactorily for the varia- 
tion of the refractive index with temperature over the whole spectral range 
from the extreme ultra-violet up to the red, and also for the course of the 
refractive index-temperature curve for a particular wavelength. 


Fused quartz is particularly interesting in these studies becavse its 
thermal coefficient of expansion is extremely small. As is well known, the 
density of quartz falls down on fusing, but thereafter, it undergoes very 
little change with temperature when the virteous silica is heated or cooled. 
Consequently, the variation of refractive index with temperature is mainly 
caused by the variation in the dispersion frequercies, the contribution cue 
to the variation in the number of dispersion centres being comparatively 
small. Ona the analogy of the result obtained with diamond, one may 
assume that the dispersion frequencies of silica also move to increasing 
wavelengths with rise of temperature. This has in fact been verified for 
the lattice frequencies (Nedungadi, 1940). Then, the contribution to dn/dt 
by this cause is positive and since the volume change is small, the negative 
contribution to dn/dt from this cause will be much smaller than the positive 
component. Consequently, one would expect that vitreous silica would 
exhibit a resultant positive dn/dt. This is in fact found to be the case over 
the whole spectrum (Sosman, 1927). 


In crystalline quartz, on the other hand, the thermal expansion is 
appreciable, and there is the possibility of the negative component of dn/dt 
being the larger one. In fact, both the ordinary and extraordinary indices 
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of a-quartz diminish with rise of temperature in the visible and near ultra- 
violet (Sosman, Joc. cit.). In the far ultra-violet they increase with rise of 
temperature, since the positive component of dn/dt becomes large owing 
to the approach to the dispersion frequencies. 


2. CALCULATION OF THE THERMAi VARIATION OF REFRACTIVE 
INDEX FOR DIFFERENT WAVELENGTHS 


The dispersion of vitreous silica has been studied by Martens and 
Micheli (1904) who have fitted their data to the following formula with one 
dispersion frequency in the ultra-violet: 


n® = 1+36112 + 0-74655 A?/A? — (0- 107044)? — 0-01350A? (l) 


the wavelength being expressed in microns. The author has verified that 
this formula fits the refractive index data from 1850 A.U. to 6500A.U. It 
may be noticed that the dispersion frequency is at 1070 A.U., which is very 
close to that used for crystalline quartz by other workers (Coode-Adams, 
1928; Radhakrishnan,1947). In the work of Radhakrishnan, two frequen- 
cies at 1060 A.U. and 1190 A.U. are employed, the latter being much weaker. 
These frequencies have actually been observed for quartz by Tousey (1940) 
as feflection maxima analogous to the Reststrahlen in the infra-red. For 
fused quartz, Tousey (loc. cit.) found two bands in very nearly the same 
positions as for crystalline quartz, but very much broader. Considering 
the relative strengths of the two frequencies, the representation of their total 
effect at regions of the spectrum not close to them by a single term employ- 
ing a frequency of 1070A.U. appears justifiable. The infra-red term in 
the above formula (1) is represented by a term of the type — KA* and does 
not explicitly involve the frequency. However, there are reasons to believe 
that none of the frequencies in the infra-red are nearer than 8y. In the 
dispersion formula of quartz (Rubens and Nichols, 1897; Radhakrishnan, 
loc. cit.), there is a weak term with the dispersion wavelength in the neigh- 
bourhood of 8 and a strong one at about 214. Upto 8p, the absorption 
of fused quartz follows very closely that of crystalline quartz, except that 
the peaks are broader (Drummond, 1936). Also, the Raman spectrum of 
fused quartz exhibits a great similarity to that of crystalline quartz, the 
lines in the latter being replaced by broad bands in the former (Gross and 
Romanova, 1929; Kujumzelis, 1935). These show that the infra-red 
dispersion wavelength A, of vitreous silica is not less than 8. Usirg a 
value of A, = 8p and 12,, it has been calculated that its contribution to 
dn/dt over the range of wavelengths we are interested in is extremely small 
compated with that of the ultra-violet frequency and that the contribution 
progresssively diminishes as A, is taken further away in the infra-red. One 
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need not therefore be disturbed at the non-occurrence of A, explicitly in the 
dispersion formula. 
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Eq. (1) may be put in the form 
n® =C + a,d,? 4/(A® — A,2) — KA® (2) 
where C =- 1-36112, a,A,? = 0-74655, A, = 0-107044, K = 0-01350. 
Differentiating (2) with respect to ¢t, and putting da/dt = — ya, 
dn y art | 1 aAPzr 1 dC 


dt ~~ MQ—2%) * nO =A)? * * In de 
where x and y have the same significance as in Part I, viz., x = (1/A,) (dA,/dt) 
and y is the coefficient of cubical expansion. The only unknown quantities 
in the right-hand side of (3) are x and dC/dt, which we may denote by 8. 


Measurements of dn/dt for fused silica at a mean temperature of 60° 
have been made by Martens (1904) for a series of wavelengths from 1850 
A.U. to 6000 A.U. Using any two measurements, x and 8 can be calcu- 
lated and then the values of dn/dt for a series of wavelengths can be calcu- 
lated. In this way, the continuous curve in Fig. 1 has been plotted showing 
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Fig. 1. Variation of dn/dt with wavelength 
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how dnj/dt would vary with wavelength according to theory. The values 
used for x and B wete x = 24:77 x 10°* and B =— 10-8 x 10". The 
citcles in Fig. 1 represent the experimental data of Martens (loc. cit.). The 
triangles represent the data of Tilton and Tool (1929). It will be seen 
that the theoretical curve fits the data well over the whole range of 
wavelengths. 


3. CALCULATION OF THE COURSE OF THE REFRACTIVE 
INDEX-TEMPERATURE CURVE 


As in the case of diamond, one may attempt to compare the value of x 
for the dispersion frequency with that for other characteristic frequencies. 
As already mentioned, the Raman spectrum of fused quartz consists of 
broad bands but in general these bands are in the same positions as the sharp 
lines in the Raman spectrum of crystalline quartz. For a-quartz, the sharp 
and intense line at 465 cm.-! stands out prominently in the Raman spectrum. 
The variation of the frequency of this line with temperature has been studied 
by Nedungadi (1940) from — 180° upto 500°C. From his data, x at 60°C. 
has been calculated to be 3-2 x 10-°. It is seen that this is of the same order 
as the value of x for the dispersion frequency of fused quartz. In fact, the 
latter is 77% of the former. 


Oae may attempt an extrapolation and assume that this relation between 
the x’s of the two frequencies is true at other temperatures, i.e., that the 
proportionate change in frequency of the 1070 A.U. band of fused quartz 
is 0-77 times that of the 465cm.* line in crystalline quartz. Then, it is 
possible to calculate the variation of refractive index with temperature over 
a range of temperatures, using (3). Here, y is known at various tempera- 
tures and x can be calculated from the above assumption. For dC/dt, 
some further consideration is necessary. As mentioned in Part I, the term 
C arises from the constant contribution to refraction of frequencies in the 
very remote ultra-violet. dC/dt thus arises from the variation in the num- 
ber of centres and frequency in this term. It may be shown that the 
contribution to dn/dt due to a dispersion frequency at wavelength A, where 
A, << Ais (a,A,?/n) (x, — y/2). Since dC/dt is negative, it follows that x, is 
less than y/2 i.e., less than 0-6 x 10° at 60°C. Thus x, is less than 2% 
of the value of x for the 1070A.U. dispersion frequency. This is not un- 
expected for these remote ultra-violet frequencies arise from transitions of 
the inner electrons and are not likely to change with temperature. In 
calculating dC/dt therefore x, was put equal to zero and dC/dt was taken 
to vary as y. 
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If An is the change in refractive index for a change of temperature At, 
then from (3) one has 


| a, A,? A? AV 1 a, A,? 4 Ar, '& 

atm ~—" O? = Vv ve a? 7 iS? ov | nd a 4 
An was calculated successively betwe@h the temperatures — 150°, 0°, 150°, 
300° and 450°. AV/V was calculated from the data on thermal expansion. 
Av/v was taken to be 0-77 times the corresponding value for the Raman 
line 465cm.-', which was determined from Nedungadi’s data (Joc. cit.) 
Since the accuracy of these data was only 1 cm.~}, a smooth curve was drawn 
through them and the change in the wave-number determined from the 
curve. The change from the value at 0° to various temperatures found 
in this way were: for — 150°, ~—1-9cm-!; 150°, 2-5cm-!; 300°, 5:3 cm-; 
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Fic. 2. Refractive index-temperature curves for four wavelengths. The continuous curves 
are from theory and the circles represent the experimental data. 
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450°, 8:7cm.+ The last term was calculated with the assumption just 
mentioned. In this way, refractive index-temperature curves have been 
computed and plotted as continuous curves in Fig. 2 for four wavelengths. 
Exderimental data for these are available from the work of Rinne (1914) 
computed to absolute indices by Sosman (loc. cit.). These are shown by 
circles in Fig. 2. The theoretical curves have in every case been made to 
pass through the measured value at 18°C. It will be seen that the theory 
fits the data for higher temperatures upto 500°C. Below 0°C., however, 
theory gives systematically a lower value than experiment. The discfe- 
pancy may be due to the fact that the dispersion frequency varies less in 
proportion to the infra-red frequency than the ratio 0-77 at these tempe- 
ratures. It may be mentioned that the error may also be partly due to the 
fact that Nedungadi’s data do not include any measurements in between 
— 190° and 0°C. 

In conclusion, I wish to thank Prof. Sir C. V. Raman for the encourage- 
ment he gave me during the investigation. 


SUMMARY 


The theory outlined in Part I of these series has been applied to calcu- 
late the thermal variation of refractive index of fused silica. It is found 
that the theory can account quantitatively for the experimental values of 
duldt over the range of wavelengths from 1850 A.U. to 6000A.U. Using 
the fact that the frequencies of vitreous silica are practically the same as 
those of crystalline quartz, the variation of the refractive index with 
temperature from — 130°C. to 500°C. has been calculated and found to 
fit well with the experimental data. 
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1. INTRODUCTION 


In Parts I and II of these series (Ramachandran, 1947), a general theory 
has been developed for the variation of refractive index of solids with 
temperature and has been applied to the cases of diamond and vitreous 
silica. As shown in these papers, the temperature variation of refractive 
index (d1/dt) can be supposed to consist of two parts, one due to the change 
in d2nsity of the solid and the other due to the changes in the dispersion 
freqiencies, the former being negative and the latter positive, in general. 
Tae thermal expansion of both diamond and fused quartz is very small, 
so that dn/dt for these is positive over the whole spectrum from tke far 
ultra-violet to the red. Fluorspar, on the other hand, has an appreciable 
thermal exsansion and exhibits a negative dn/dt throughout this region. 
Tae value is about —5~x 10 at 1850A.U., which increases numerically 
to about — 12x 10 at 6000A.U. (Micheli, 1902) and then again 
diminishes to about —9x 10 at 6-5 (Liebreich, 1911). It will be 
interesting to see how far the theory can explain these variations. 


The dispersion of fluorspar has been studied by various workers, and 
a number of dispersion formule have been suggested (Jaffé, 1928). The 
following formula given in the International Critical Tables, Vol. VU, p. 12, 
is said to fit the experimental data from 0-185, to 9-4u: 


0-00612093 5099-15 

nt= 6-09104 + s— 9.008884 * X¥— 1258-47 @) 
This makes use of an ultra-violet frequency at A942-6A and an infra-red 
one at A35-48y4. The ultra-violet frequencies of fluorspar have recently 
been investigated by Tousey (1936) using a technique analogous to that used 
for the Reststrahlen in the infra-red. It is found that there are two 
absorption bands, a narrow one with a half-width of roughly 100AU. 
with its centre at A1115A.U. and a stronger and more diffuse band at 
wave-lzngths shorter than 1060A.U. and having its maximum at about 
909 A.U. or slightly shorter wave-lengths. In view of Tousey’s observa: 
tions, it will be seen that the above formula (1) can at best only te an 
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approximation. In fact, it is found that below 1850 A.U., the formula 
progressively departs more and more from the experimental values of 
refractive index as one goes into the ultra-violet. The transmission limit 
of fluorspar is about 1250 A.U. and refractive index data down to 1311 A.U. 
are available from the work of Handke reported by Martens (1923) in the 
Landolt-Bornstein Tables. The author has been able to fit these data as 
well as the data in the visible by a single formula of the form 


nia C+ BAA(e — i,3) (2) 


with two ultra-violet frequencies at A1115A.U. and 820A.U. and an 
infra-red frequency at A35-48u4. In fact, a formula of the type 


(nt— 1) = ZA,A4/(A*® — 2,2) (3) 


can also be derived having four dispersion frequencies at A450A, 888 A, 
1115A and 35°48. Details of these will be given in later sections. 


The question now arises as to which of these formule has to be used 
for calculating dn/dt. It will be evident that as long as a certain formula 
fits the data for refractive index in the range of wave-lengths in which 
one is interested, then the theory developed in PartI can be used with that 
formula for that same range. Consequently, formula (1) has been used for 
the purpose, which has the advantage of simplicity. It is found that it 
explains quite well the variations in dn/dt described earlier. Later, formule 
(2) and (3) have been applied, and they lead to certain interesting results. 
It is found that the extreme ultra-violet frequency at A450 A does not vary 
with temperature, the one at A888A varies only slightly, while the one 
at A1115A varies much more with temperature. 


2. DISPERSION FORMULZ FOR FLUORSPAR 
Formula (1) mentioned in the previous section can be put in the form: 


4 0-68898 2, 4-05186 A? 
PI: nt= 1-33020 + ss="9-008884T NF 1258-47 4 


Two new formule have been developed by the author of the types 


(2) and (3). They are: m ni a “nee 
0-61667 16404 -051 
PIL: n= 1-25910 + y5—9.9067240 X= 0-012432 7 F— 1258-47 ©) 











045848 a? 0-4282442 . -0-15297 a8 
PII: n*— 1 = y3—9-9020250 + 18—0-0078925 * 1®— 0-012432 
+ 4:05186 x - 
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and are meant to represent the refractive indices m with respect to vacuum 
at A, the wave-léngth in vacuum. In Table I, the difference between the 
values calculated from the above formule ard the experimental data are 
given as the number of units in the fourth decimal place. (In Marten’s Table 
quoted earlier, the refractive indices have been given only to 4 decimal 
places.) It may be mentioned that F 1 refers only to refractive indices in air. 
Making correction for the refractive index of air, the errors are shown in 
brackets under FI. It will be seen that the error becomes larger and 
larger as one progresses into the ultra-violet beyond 0°185y. FIL and 
FIII on the other hand fit very well down to 0:16 beyond which both 
show a peculiar type of error. The difference between experiment and 
theory is first positive, becomes larger artd larger, reaching a maximum at 
about 0-145, then diminishes and becomes negative. The author fouud 
that no manipulation with the frequencies below A 1115 A.U. can improve 
the fit in this region. In fact, when the experimental dispersion curve is 
plotted, an anomalous rise in refractive index is found down to about 
0-145 and then an anomalous fall is found. A similar variation is found 
in Tousey’s reflection curve also (loc. cit.). All these suggest that there 
is an extremely weak absorption band with its centre at about 0-134, (at 
which waye-length the error passes from positive to negative values), which 
is the cause for the anomalous dispersion in its neighbourhood. 




















TABLE | 
nexp—"th 

: Avac *tpac ©XP- 

in microns FI | FII FIII 
0- 76845 1-4313 0 0 
0-58947 1 +4342 +3(-) 0 0 
0-394624 1 -4427 0 0 
0-308314 1 +4530 0 0 
0-257392 1+4653 +4 (-1) 0 0 
0-219527 1-4821 0 0 
0-199057 1 -4969 +8 (+2) 0 0 
0-185477 1-6107 +10 0 0 
0-176420 1 +5227 +14 0 0 
0+164455 1-5438 0 0 
0-15692 1-5615 —J -1 
0-16234 1-5751 +46 +1 +1 
0-14817 1-586 +2 +2 
0-14363 1.6089 +88 +6 +5 
0-13940 1-6305 +4 +3 
0°13536 1-6565 +1 +1 
0-13192 1-6844 +237 -3 -2 
0-13111 1-6921 
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3. CALCULATION OF dn/dt FOR VARIOUS. WAVE-LENGTHS 


A. Using FI—As was mentioned earlier, formula FI represents the 
refractive index sufficiently well for wave-lengths larger than 1850A.U. 
One may therefore calculate dn/dt from it inthe same manner as was done 
for fused quartz. If one writes 


nt=C + Za, A*, Ah A®— A,%), (7) 
“an dn 8B ( “ A,? »# 
y(n?—C) , 1 3 Or Xr 
a «wa. +; “ — 4,7)? (8) 


where B =(dC/dt), x, =(1/A,) (da,/dt) and y is the coefficient of cubical 
expansion. The value oi xs for the infra-red frequency is available from the 
unpublished work of D.C. Press carried out in this laboratory on the 
temperature variation of the principal Raman frequency at 320cm- It is 
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Fic. 1. Curves showing the variation with wavelength of the contributions to dn/dt due 


to volume and frequency changes. The abscisse have been reduced five-fold after 1 « to facilitate 


the whole range from 0-2 4 to 7 u being shown in the same diagram. The scale of the ordi- 
nates for both the curves is the same, but the origins have been shifted, 
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found that the contribution to dn/dt due to the change in the infra-red 
frequency is negligible for the visible and ultra-violet regions of the 
spectrum, being less than 1% at 0-8u and less at shorter wave-lengths. 


The calculations were done for 60°C. for which temperature the dn/dt 
data are all available. y was obtained from the unpublished work of 
D C. Press in this laboratory. There are only two unknown quantities, 
x, and 8, which were calculated to have the values 


X= 33-1 x 10°* and B= — 37-6 x 10 (9) 


As mentioned earlier, dn/dt can be supposed to be composed of two 
parts, due respectively to the change in density and change in dispersion 
frequency, which we may denote by P and X, following the notation in 
Part I. Thus, 

on on da 

pam eK (10) 
Now, in the expression (8) above, the term 8/2 n is the contribution due to 
the alteration of the constant C, which is negative. Following the reason- 
ing advanced in Part II with respect to vitreous silica, one may conclude 
that this part is wholly due to the change in volume. In other words, it is 
assumed that the extreme ultra-violet frequency which gives rise to a 
constant term in dispersion does not vary with temperature, which is 
a reasonable assumption. Thus, 


P=— y (n?—C)/2n+ B, X=(1/n) Za, d,*x,[(A*— A,*)* (11) 


In Fig. 1, these two parts are shown separately by the curves (a) and (6). 
The calculations for drawing the curves were made with the values for 
y, x, and B given above. For the infra-red, x, was taken from Press’s work, 
which at 60°C. has the value 159x 10-*. The curves show interesting 
variations of P and X. The contribution due to volume change follows 
a course similar to the refractive index, diminishing (algebraically) as one 
approaches the ultra-violet frequency and increasing towards longer wave- 
lengths in the infra-red. The other contribution, however, increases on 
approaching a dispersion frequency, the increase taking place on either 
side of it. This is because it is inversely proportional to (A?— A,?),? and for 
the same reason, the increase in X preponderates over the diminution in 
P as one progresses into the ultra-violet. Thus, dn/dt increases on approach. 
ing both the ultra-violet and the infra-red dispersion bands. In Fig. 2 the 
course of the (dn/dt)— A curve is shown, which will illustrate this fact. 
The continuous curve has been calculated from theory, while the points 
are the experimental data of Micheli (1902) and Liebreich (1911), the former 
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being shown by circles and the latter by triangles. It will be seen that the 
agreement With Micheli’s data is remarkably good. Liebreich’s data exhibit 
an appreciable algebraic increase in dn/dt from 0-8 to 2u and a relatively 
smaller increase thereafter than indicated by theory. In this region, theory 
ives systematically lower values, although it agrees with the data at longer 
wave-lengths. The author is unable to conceive of any reason as to why 


gnjdt should quickly increase (algebraically) upto 24. In fact, as will be. 


seen later in this section, the use of other dispersion formule like FII and 


FIII does not alter the theoretical curve sensibly and that the discrepancy, 


is still present. It would appear that the discrepancy between theory and 
experiment in this region must be attributed to experimental errors, which 
are likely to be large since dn/dt is only obtained indirectly from the difference 
in the refractive indices at two temperatures and also since the latter 
measurements themselves are less accurate than in the visible because no 
photographic method can be used. 
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——> Wavelength in microns 


Fig. 2. Variation of dn/dt with wavelength. The circles are the experimental data of 
Micheli and the triangles those of Liebreich. The scale has beem compressed five-fold above 
1:0y. 

B. Using F II and F IlIl—We lave thus seen how a relatively approxi- 
mate formula like FI is sufficient to explain the variations in dn/dt. We 
shall now apply the new formula FII, which is also of the form (7), so 
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that dn/dt is given by (8). Now, there are two ultra-violet frequencies at 
\1115.A.U. and 820A.U. For convenience, we shall assume that both of 
them vary with temperature to the same extent in proportion to their 
frequency, in other words that both have the same value of x. Then, it is 
found that B and x Have the values: 

x =33-0x 10°, B=—37-6x 10° (12) 
The resulting formula fits the data for dn/dt very well. In Table II below, 
the calculated values of dn/dt for a few wave-lengths from FI, FII and 


FIII are given together with the experimental values. It will be seen that 
there is little to choose between the three formule. 


However, FIII is. interesting because it is of the standard type given 
by the Ketteler-Helmholtz theory with the constant C= 1 in (7). Conse- 
quently, dn/dt does not involve an arbitrary constant B, but depends only 
on the values of x, for the dispersion frequencies. Denoting the four 
frequencies at AA450A.U., 888A.U., 1115A.U. and 35-47 respectively 
by the subscripts 1, 2, 3 and 4, one finds that the x,’s have the values: 


X= 0, a 11:3x 10, Xs= 70:2 x 10-* (13) 


x, Was taken to have the value 159 x 10~* obtained from Press’s measure- 
ments (loc. cit.). 











TABLE II 
dn/{dt calculated from 
Wavelength in dn/dt from 

microns experiment 

FI FII FIII 
0+ 199057 — 6°07 | — 5-92 | — 5-93 — 5-89 
0+257392 — 9°26 | — 9-19 | — 9-39 — 9°24 
0-308314 —-10°-45 | —10°36 | —10-44 — 10°32 
0-68947 ~11-84 | —11-91 —11-78 — 11-93 
1-3 — 12-23 —12-16 | —11-91 —11-21 
3°16 —11-63. | —11-55 —11-28 — 9°83 
6-5 - 9-15 | — 9-05 | — 8-78 — 8-84 

















It will be seen the proportionate rate of change of the ultra-violet frequencies 
becomes less and less as one proceeds further into the ultra-violet. Thus 
the remote dispersion frequency at A450 A.U. is unaffected by temperature 
while between the two at 1888 A and A1115A., the latter increases much 
more than the former. It may be mentioned that the actual numerical 
values given in (13) can be altered slightly and yet fit the experimental] data. 
However, it will be found that always xs> x2> x, the last one being very 
small. 
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3. DISCUSSION OF THE RESULTS 


The result just mentioned is what one would expect from other 
considerations. Of the various ultra-violet frequencies, those having higher 
frequencies should arise from the deeper levels of the atoms, while those 
with lower frequencies arise from transitions between the levels of the 
outer electrons. Since the outer electrons take part principally in the binding 
of the crystal, a change in temperature would be expected to affect the 
levels of the outer electrons much more than those of the inner ones. 


Another interesting fact emerges when one compares the value of X 
for the ultra-violet frequencies with that for the characteristic infra-red 
frequency of the lattice. It was shown in Parts I and II that these two are 
of the same order in diamond, and that in fused quartz x of the dispersion 
frequency is of the same order as x of the intense Raman line at 465 cm 
of a-quartz. For fluorspar, although the temperature variation of the infra- 
red active lattice frequency has not been studied, data on the Raman 
frequency at 320cm.-! (which lies close to the former) are available (Press, 
loc. cit.). If one takes FI or FII, then x for the ultra-violet frequencies 
is about 33 x 10-* which is only about 1/5 of x for the infra-red, viz., 
1599x 10-*. If we denote by xy and xp the values of x for the ultra-violet 
and infra-red frequencies, then xy/xp is 78% for diamond and 77% for 
fused quartz while it is only 21% for fluorspar. Even if one takes FIII and 
considers ys, the variation of the frequency at A1115A.U., then also, 
yalXe is only 45%. Thus, the ultra-violet frequencies vary much less 
relatively in fluorspar than in diamond and fused quartz when compared 
with the variation of the lattice frequency of the respective crystals. The 
interpretation of this difference is not far to seek. Both diamond and 
quartz are typical covalent solids, in which the binding is éssentially 
produced by the sharing of the outer electrons. One would therefore 
expect the changes in the binding energy and changes in the electronic 
states of the outer electrons to correspond closely in these solids. In 
fluorspar, the binding is not so highly covalent as in diamond or quartz, 
but is partly electrovalent. It is therefore not surprising that the ultra-violet 
electronic transitions are relatively less affected by temperature than the 
infra-red lattice frequency. 


One may also consider what happens to xy at elevated temperatures. 
For this, one can make use of the work of Reed (1898) who has determined 
dn/dt for sodium light for various temperatures above the room temperature. 
y Can be obtained from the data of Press mentioned earlier. It is assumed 
that the component of dn/dt produced by change of volume varies.in the 
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same proportion as y. This means that 8 is taken to vary as y, which is true 
since f is due only to a change in volume, the extreme ultra-violet frequency 
not altering with temperature. The important values used in the calcula. 
tion, as well as the values of xy (calculated with FI) and xg (calculated 
from Press’s data) are shown in Table III. 

















TABLE II1 
Temperature °C | dn/dtx10®| yx10® XyX10° | x_,x10° 
| 
60 -11-9 | 63-24 33-1 159 
100 -12-4 | 64-17 32-8 | 181 
150 -18-1 | 65-58 32-7 | 110 
200 -13-7 | 67-50 33-2 | 04 
250 | ~14+3 | 10°35 34-4 90 
| 





It will be noticed that yy is practically a constant upto 200°C. after which 
it rises slightly. yg on the other hand shows an appreciable decline with 
rise of tempearture. In both diamond and fused quartz, xy increases with 
temperature. Although the precise reason for the difference in the case of 
fluorspar is not clear, one may reasonably attribute it to the differences in 
the nature of binding in the crystal. 


I am very grateful to Prof. Sir C. V. Raman for the kind and encourag- 
ing interest that he took in this investigation. I am also thankful to 
Mr. D. C. Press for having kindly made available to me the unpublished 
results of his studies on fluorspar. 


SUMMARY 


Making use of the dispersion formula given in the International 
Critical Tables, with absorption bands at AA0-0942 and 35-48 y, and using 
the theory developed in Part I, the course of the variation of dn/dt with 
wave-length has been satisfactorily explained from 0-185 to 6-5y. The 
dispersion formula however fails below 0-185, and consequently two new 
formule have been developed holding down to 1300A.U., one of which 
uses dispersion frequencies at AA 0-082, 0-1115 and 35:48. together with 
a non-unity constant in the expression for m?, while the other is of the 
Ketteler-Helmholtz type with four frequencies at AA 0-045, 0-088, 0-1115 
and 35:48. All the formule explain the variations of dn/dt, in particular, 
as to why it increases algebraically as one proceeds both into the ultra- 
violet and the infra-red. The interesting fact emerges that the extreme ultra- 
violet frequency at A0:0454 does not vary with temperature, while the 
proportionate variation of the one at 0-0888, is much less than that of the: 
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one at 0°1115u. This result becomes intelligible when one remembers 
that the deeper levels in the crystal would be less affected by temperaiure 
than the low-lying ones. Even the last one, although it is of the same 
otder, is only one half of the proportionate variation of the infra-1ed 
lattice frequency. It is also found that the rate of variation of tke vltra- 
viol2t frequencies is constant at different temperatures. In these respects, 
fliorspar differs from diamond and vitreous silica. The difference has been 
attributed to the fact that the binding in fluorspar is not purely covalent, 
but is also partly electrovalent. 


REFERENCES 
1. Jaffé, G. .. Handb. d. Exp. Phys., 1928, 19, 176. 
2. Liebreich, E. .. Verh. Deutsch. Phys. Ges., 1911, 9, 1. 
3. Martens, PF. F. .. Landolt Bérnstein Tables V Edn., Hauptwerke, 1923, 
Vol. II, 911, 
4. Micheli, F.J. .. Ann, d. Phys., 1902, 7, 772. 
5. Ramachandran, G. N. .. Proc. Ind. Acad. Sci., A, 1947, this issue, 
6. Reed .. Ann. d. Phys. u. Chem., 1898, 65, 734. 
7. Tousey, R. Phys. Rev., 1936, 50, 1057. 











94-47 Printed at Tne Bangalore Press. Bangaiore City by G. Srisivasa Rao, Superintendent 
and Published by The Indian Agademy of S Bangalore 








Description of Plate X, Vol. XXV, A, 1947 


1. Spectrum (a) is tek2n on E, Hilger quartz spectrograph using 


‘Kcedak Superpanchro press P-1200’ plate. 


2, Sp:ctrum (5) is takzn on M:d‘um Hilger quartz spzctrograph usirg 
‘Kodck Process Regular B-20’ plate. 


ERRATA 


Vol. XXV, A, No. 2, February 1947 
First liae in foot rote, for {[o,*, I7,§, IT,5 "IT ,3q] o,*} 
read {[o,%, 7,4, 7, + *IT4] 9,*} 
Fourth line, for 1, (°2}) read 1, (*25) 
23rd line, for O} (*2,) read O¢ (*2-) 
28ih lire, for 51583 read 51528 





